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Theorie des Gasdrucks gegen eine bewegte Wand 
Be Von @. Schweikert in Dingelstädt (Eichsfeld) | 


Wand entwickelt, die der Änderung des Impulses und der Geschwindigkeitskomponente der stoßenden Mole- 

 küle in Abhängigkeit von der Geschwindigkeit der Wand Rechnung trägt, aber gegebene Dichte und Tempe- 
ratur des Gases unmittelbar an der bewegten Wand voraussetzt. Es wird eine wichtige Formel gewonnen, 
die zeigt, daß der Gasdruck für ein ideales Gas mit wachsender Wandgeschwindigkeit asymptolisch gegen 
Null geht, während für die realen Gase der Gasdrusk bereits bei endlichen Werten der Wandgeschwindigkeit 
— je nach .der Größe der Abweichung vom idealen Gaszustand — auf Null herabsinkt. 


E | | 
E- Basing on the kinetic gas theory, a theory of the gas pressure against a moved wall is developed that 
takes account of the change of ihe impulse and of the component of velocity of the impact molecules but presumes 
; the knowledge of ihe density and ihe temperature of ihe gas close by the moved wall. An important formula 

is obtained proving that while the Pressure in- case of an ideal gas converges to zero with increasing velocity 
öf the wall ihe pressure of a real gas is already declining io zero at finite values of the velocity of Ihe wall 
according to Ihe amount of deviation from the state of Ihe ideal’ gas. 


Mexona us KUHeTnyecKoü TeOopuHH Tas0B PaspabaTblBaerca Teopua MaBleHng Tasa Ha 
YAapAIOINUX MONeKYI B 3&BUCHMOCTH OT CKOPOCTH CTeHEl, & OTHOCHTENBHO INIOTHOCTH .U 
TeMHEePaTyPpEI Ta83 IIPeADONOTAETCH, YTO OHH MAHEI HEIIOCPeNCTBEHHO Yy ABHIKYINeÄCH CTEHB, 
. BurBonuTea cyıNecTBeHHaa hopMyAa, U3 KOTOPoü BHAHO, UTO NABIIEHHe HNEAJIBHOTO TA3a ITPu 

BOSPACTalmINea CKOPOCTH CTEHEI ACHMTOTHYECKH CTPeMHTCH K HYIII, TOTNA KAE B CAIyyae 
Bu: HeNNeaI5BHOTO TA33 NABIeHHE TANAET NO HYJIA YıKe IIPH KOHEYHEIX CKOPOCTAX, COOTBETCTBEHHO 
CTeIICHH OTKJIOHCHHA OT HNEANIBHOTO COCTOAHHA TA3A, 


Einleitung 


2.2, Bereits vor 25 Jahren habe ich zum ersten Male darauf, hingewiesen 1), daß aus der kine- 
- tischen Theorie der Gase eine Änderung des Gasdrucks gegen eine bewegte Wand gegenüber dem 


-  Gasdruck auf eine ruhende Wand folgt, und habe damals eine einfache Formel für die Abhängig- 
{ keit des Gasdrucks von der Geschwindigkeit der Wand unter vereinfachenden Voraussetzungen 
. abgeleitet. Die damals entwickelte Theorie trug zunächst der Anderung des Impulses Rechnung, 


der von den stoßenden Molekülen auf die bewegte Wand ausgeübt wird. 


Seitdem sind noch einige Abhandlungen?) erschienen, die gleichfalls nur die Änderung des 
Impulses bei gleicher mittlerer Molekulargeschwindigkeitskomponente ohne Annahme eines Ver- 
teilungsgesetzes und nur den idealen Gaszustand berücksichtigen, oderaberauf hydrodynamischer 
Grundlage die Änderung der Dichte an der Wand infolge der Strömungsvorgänge in den Gasen 

zu berechnen versuchen. | 


Es sind nun drei grundsätzlich verschiedene Umstände, welche die Änderung des Gasdrucks 
auf eine Wand infolge ihrer Bewegung bedingen und zwar ist es: 1. die Änderung des Impulses, 
der von den stoßenden Molekülen auf die Wand übertragen wird, infolge der Bewegung der 
Wand; 2. die Änderung der Geschwindigkeitskomponente der stoßenden Moleküle nach der 
| Reflexion an der bewegten Wand, — diese Änderung der molekularen Geschwindigkeitskom- 

ponente senkrecht zur Wand infolge ihrer Bewegung ist zugleich die mechanische Ursache für 
“ die Änderung der Temperatur des Gases bei adiabatischer Zustandsänderung durch Arbeits- 
leistung gegen einen Kolben —, und 3. die Anderung der Gasdichte an der bewegten Wand 


infolge der Änderung des Volumens, das der freien Bewegung der Gasmoleküle zur Verfügung. 


steht. In der folgenden Theorie zur Ableitung des Gasdrucks gegen eine bewegte Wand wird 
die dritte der vorstehend genannten Ursachen zunächst außer acht gelassen. Ferner berück- 
sichtigen wir von der zweiten der obigen Ursachen für die Änderung der Geschwindigkeits- 
komponente der Moleküle durch den Stoß gegen die bewegte Wand zunächst nur die Wirkung 
dieses Umstandes auf die Änderung des Gasdrucks, nicht hingegen die Wirkung auf die Anderung 


| 1) Vgl. Innere Ballistik. Leipzig 1923, 8. 58. - 2 
2) Da mir keine Literatur zurzeit zur Verfügung steht, kann ich diese Arbeiten leider nicht näher bezeichnen. 


19 


e: j Band ZI 5 Oktober 1949 i N Heft 10 


Auf der Grundlage der kinetischen Theorie der Gase wird eine Theorie des Qasdruckes gegen eine bewegte 
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der Temperatur des Gases. Wir betrachten also in unserer nachstehenden Theorie die Dichte 
und Temperatur des Gases unmittelbar an der bewegten Wand als vorgeg bene 
Größen. |  E 


$1. Ableitung des Druckes, der von den Molekülen REN 
einer gegebenen Geschwindigkeitskomponente (x) auf eine bewegte Wand ausgeübt wird 

Wir gehen in unserer gaskinetischen Theorie davon aus, daß der von den stoßenden Mole- 
külen auf die bewegte Wand übertragene Impuls nach den Gesetzen des zentralen, elastischen 
Stoßes eine Änderung in Abhängigkeit von der Geschwindigkeit (v) der Wand erfährt ‚und 
zwar ändert sich nach diesem Gesetz die Geschwindigkeitskomponente (x) des stoßenden 
Moleküls in der zur Wand senkrechten Richtung. Die Geschwindigkeitsänderung (dv) der Wand 
infolge des Impulses eines einzelnen mit der Geschwindigkeitskomponente (@) senkrecht zur 
Wand stoßenden Moleküls beträgt bekanntlich®): 


d 


v— am (e — v) 

 M-+m ; 

wobei m die Masse des Moleküls und M die Masse der Wand ist. Da letztere stets sehr groß 
gegen erstere ist, kann diese Gleichung einfacher geschrieben werden: 


wel... ..- 02.00 


Die durch den Stoß auf die Wand ausgeübte Kraft beträgt: 


wenn o die Oberfläche der Wand, p den Druck und di das (unendlich) kleine Zeitelement be- 
zeichnet, in dem im Mittel ein Molekül mit der Geschwindigkeitskomponente x gegen die Wand 
stößt, oder anders ausgedrückt: dt ist das mittlere Zeitelement zwischen zwei aufeinander 
folgenden Molekülstößen mit der Komponente x gegen die mit der Geschwindigkeit v bewegte 
Wand. 

Die wesentliche Aufgabe unserer Theorie besteht nun darin, dieses mittlere Zeitelement dt 
in Abhängigkeit von allen dasselbe beeinflussenden Momenten der Molekülbewegung zu be- 
stimmen, also insbesondere unter dem Einfluß der mittleren freien Weglänge (Z), des mittleren 
Radius (o) ‚der mittleren Stoßdauer (r) der Moleküle, sowie unter dem Einfluß der Geschwindig- 
keit (v) der Wand. 

Wir berechnen zunächst, wie auch früher, das gleichfalls praktisch unendlich kleine Zeit- 
intervall dt,, in dem im Mittel die Anzahl n, der in einer Elementarschicht von der Höhe 2, 
vorhandenen Moleküle mit der Geschwindigkeitskomponente x in Richtung der mit der Ge- 
schwindigkeit v bewegten Wand gegen diese stoßen. Entsprechend diesen Festsetzungen ist 


DAmEir Rn ar A (3). 


Wir betrachten nun ein bestimmtes Molekül der Geschwindigkeitskomponente x in der 
Elementarschicht }, an der bewegten Wand. Sein Abstand von der Wand im gegebenen Augen- 


: 1 \ ar A 
blick betrage 7 l,, wobei der Wert 21 sei. Die Zeit, die vergeht, bis es gegen die Wand stößt, 
sei dt,. Dann gilt die Beziehung: 


DE ee Ge ER ve a Pu Ne N a 


dh: EB ER ET (4), 


da die Wand sich in der Zeit di, um die Strecke »- dt, fortbewegt. Der Faktor k hat dabei 
den Wert): 


PEN LACH. 
*) a.2.0.8.6und G.Schweikert, Zur Theorie der Zustandsgleichung der Gase. Z. Physik 1934. 
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- Er trägt der mitileren Ver küran (eo) und der mittleren Stoßdauer (7), sowie der hllden 
freien Weglänge (Z) der Moleküle Rechnung. Der Faktor k ERaIz erisiert is die ARweRhans 
des Gases vom idealen Zustand, für den k=1 ist. 


Nach dem Stoß muß sich das betrachtete Molekül bis zum Abstand (1 — =. 1, zurück- 


| bewegen, um in die der Lage vor dem Stoß entsprechende relative Lage zur Wand nach dem 
Stoß zu kommen. Die hierzu benötigte Zeit sei di,. Es gilt dann: 


ern RENTE Br 


Fa 2-— 20 


da die Geschwindigkeitskomponente naeh. dem Stoß den Wert &—2v angenommen ir infolge 


der Bewegung der Wand, und da sich außerdem der zu durchlaufende Weg um den Betrag ®- dt, 
verkleinert, wenn die Lage des Moleküls relativ zur Wand während der Dauer der Wandbewegung 
unverändert bleiben soll. Aus den beiden Gleichungen (4) und (5) läßt sich das le Zeit- 


element dt, in. Be Weise berechnen. Zunächst folgt: (e —kv)d, =— kl, und | 


(2 heldy="—I.r: > Also ist: 


TR Re 
ee Fr su E 


Durch einfache ‚Umformung erhält man daraus: 


en BE 
@_dverka—h 


Aus diesem Wert di, erhält man schließlich durch Summierung über alle möglichen r-Werte 


h di, =k-): 


‚ das gesuchte Zeitelement at,, in welchem die n, Moleküle der Elementarschicht 2,, deren Ge- 
schwindigkeitskomponente zwischen x und x +dz liegt, gegen die bewegte Wand stoßen. 


Da wir gleichmäßige Verteilung der Moleküle über den Raum voraussetzen, so läßt sich der 


Faktor 1/r in Gl. = durch folgenden Ausdruck darstellen: 
N v‚e1,233. EISEN RER nn (7), 


E 
wobei n, eine gegen die Einheit sehr (unendlich) große Zahl ist. So folgt weiterhin: 


x —k- 0 (1—h) v. S. 


di, = 1 er) 
Es ist aber: 


. Das mittlere Zeitintervall zwischen zwei HE Molekülstößen mit der 
Geschwindigkeitskomponente x beträgt somit: 


Für den idealen Gaszustand ist k = 1, und daher reduziert sich vorstehender Ausdruck 
auf die einfachere Form: 


_ Die durch den Stoß von n, Molekülen der Elementarschicht Z, mit der Geschwindigkeits- 
komponente x auf die bewegte Wand ausgeübte Kraft sei entsprechend Gl. (2) mit K, bezeichnet. 
Dann ergibt sich: 


- u. ($) ae "m 20 a 
KR=M BEE [a —2vr +k2 —k) vj=0'M: 
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Daraus folgt unter Berücksichtigung von Gl. (3): 
vl 202 +k2— MN. 
x 


N ; . 

Darin ist noch 5 = Er , während 7? = LAPe, daher = —-:P. Ferner ergibt das Verhältnis 
N, N RER £ 

(m) von Molekülmasse und Volumen .des Elementarwürfels die Massendichte : des Gases 


3 
mM 


5 kan . N, ist die Zahl der Moleküle mit der Geschwindigkeitskomponente £ im Gas- 
Pi 


volumen V, und N die gesamte Zahl der Moleküle überhaupt. Es läßt sich daher auch schreiben 


ee cn 


Dieses ist also der Ausdruck für den Druckanteil p, (Partialdruck), den die Moleküle mit der 
Geschwindigkeitskomponente & bei ihrem Stoß gegen die mit der Geschwindigkeit v» bewegte 
Wand ausüben. 


“ $2. Einführung des Max wellschen Verteilungsgesetzes 


‘Wenn aus der Gl. (10) der Ausdruck für den Gasdruck auf die mit der Geschwindigkeit v 
bewegte Wand gewonnen werden soll, so muß das Gesetz bekannt sein, nach dem die Geschwin- 
‚digReitskomponente x über die Anzahl N, der Moleküle verteilt ist. Da die Grundlagen für die 


- Aufstellung eines genaueren Verteilungsgesetzes in dem hier betrachteten Fall noch fehlen, so 


wählen wir in gebräuchlicher Weise das Max wellsche Verteilungsgesetz, beachten aber, daß 
seine Anwendung allgemein einer Idealisierung des wirklichen Gaszustandes entspricht, da die 
"Moleküle durch den Stoß gegen die bewegte Wand einseitig ihre Geschwindigkeitskomponente (x) 
ändern, wodurch eine zunehmende Abweichung vom Maxwellschen Verteilungsgesetz mit 
wachsender Geschwindigkeit (v) der Wand bedingt sein wird. 


Nach dem Max wellschen Verteilungsgesetz ist nun die Anzahl N, der Moleküle mit 
der Geschwindigkeitskomponente & bekanntlich durch folgenden Ausdruck gegeben: 


Führen wir Gl. (11) in Gl. (10) ein und setzen zugleich zur Abkürzung den Wert der relativen 
Geschwindigkeitsgrößen 


wobei & den wahrscheinlichsten Wert der Molekulargeschwindigkeit bezeichnet, so gewinnen wir 
das Druckdifferential 


Der Faktor von v/ ist hierin eine Funktion von k, er ist also durch die Abweichung des 
Gases vom idealen Zustand in seiner Größe bestimmt. 


Den Wert von k gemäß Gl. (4*) formen wir folgendermaßen um: 


Dann wird: k (2 — k) = (149) (1— 9) = 1— 92. Da in dem Ausdruck für # die Stoßdauer r 
der Moleküle offenbar noch eine Funktion der Molekulargeschwindigkeit ist, so müßte für 
eine strenge Auswertung der Gl. (13) r als Funktion von & gegeben sein. Wir können uns jedoch 
dadurch helfen, daß wir entweder r als einen konstanten Mittelwert ansehen, der nur eine Funk- 
tion der Temperatur des Gases ist, oder aber noch einfacher die Größe ® selbst als einen kon- 


'anten Mittelwert betrachten, indem wir statt & die mittlere Geschwindigkeitskomponente a 
 Substituieren und den Ausdruck a  rals eine gegebene Temperaturfunktion des Gases betrachten, 


also statt Gl. (10%) den Ausdruck setzen: | 9 — ——= |. Dementsprechend ergeben sich 


ur 


zwei ein wenig verschiedene Wege zur Auswertung der Gl. (13). Nr RE 


E ee Wir setzen weiterhin:| _ € n (ll: Rn | 5 Bi 
= ER EEE Ver a ar. {97 &T r 3 SS 2 
R L+eo: L4+eo co oa " Ei ER | re 
er, Hierin ist 7, (6) eine Korrekturgröße, die allein von der Dichte (6) des Gases blänt. 
: Für nicht zu hohe Gasdichten ist sie sehr nahe gleich dem Verhältnis 7/6 des Eigenvolumens der a 
Moleküle zum Volumen des Gases und geht mit abnehmender mittlerer freier Weglänge (Z) Eu 


gegen den Grenzwert Eins (2-01). Für hohe Gasdichten wird. aber n, eine kompli- 
ziertere Funktion der Gasdichte. | 


En 


Y 


Die Größe o = — ist hingegen eine reine Temperaturfunktion (7), die mit zunehmender 


Temperatur stets abnimmt. Unterhalb des Boyle-Punktes (75) ist ihr Wert größer als Eins, N 
während er oberhalb des Boyle-Punktes kleiner als Eins wird: ; > 
T<T;; o>1 und. a ee E Er 
Nunmehr läßt sich der Wert von © folgendermaßen schreiben: | u 

ge En ER . (16); 2 
und somit wird ' | DR | a 
1—-®=1—(02— 1 =1l—- 2 — (m'o'2)? +22 0:2. 3 


BB EZ a u han zain 030 2.0, 


Damit nimmt Gl. (13) folgende Form an: 


r h .n? j \ $ v2 
dp ’ = do) lo) tHl—moilerdz (IM. 3 
Die den zweiten obigen Festsetzungen, die sich auf die Bedeutung von ® beziehen, ent- | E 


sprechende Gleichung für dp, geht aus der vorstehenden Gleichung dadurch hervor, daß man: 

| o=(0 und |7, durch #2] ersetzt: a 
De 102 

dp, = — '— 

nee. 

Da nun von der gesamten Zahl (N) der Moleküle nur diejenigen gegen die Wand stoßen 

können, deren Geschwindigkeitskomponente & ‚größer als die Geschwindigkeit v der Wand 


(2> v,) ist, so ist der vorstehende Ausdruck über dz von v, bis oo zu integrieren, wodurch sich 
der gesuchte ‘Wert (p) des Gasdrucks gegen die bewegte Wand ergibt: 


[2 2u2 +1) dlrerdz ....... (18). 


ö:02 2 er BEI | | N 2 y 
P= —'—' za. e ?d2— 20° Zee ”dz- (1— 9) u- e z .dz , (19). 
. kg Ya v v, ® 
2 re T r 


$ 3. Auswertung der auftretenden Integrale 
Entsprechend der vorstehenden Gl. (19) sind also folgende drei Teilintegrale auszuwerten: 


© 5 1 x e ee 
I; [e;era:= 43er] +, [era r + io). 
r 
172 x Y 


Dieses Integral ist leicht zu lösen durch partielle Integration, indem man setzt: 2 e*”"d2= du; 
z=uu- . e*;dv=dz. Die Funktion © (v,) ist darin in bekannter Weise das Gauß- 


sche Fehlerintegral. 


er . £ Br 
g Be 


+ ; 
%y | 
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Die Summe aus ae 3 ER entsprechend dem Klammerausdruck in n Gl. ee Be 
zeichnen, wir mit En " B{r). Bio) ist also die Korrekturfunktion, mit der der 


19°] auf eine de Wand ae werden muß, went man den Gasdruck Pr, 


rm 3 Bio 

g 
die ınit der relativen Geschwindigkeit % -- bewegte Wand erhalten will unter sonst BIER 
Zustand des Gases an der Wand, also bei gleicher Dichte ö und gleicher Temperatur bzw. glei- ER 
chem Wert (a) der wahrscheinlichsten Molekulargeschwindigkeit. 7’ 


Wir haben damit folgende endgültigen Ausdrücke für den Gasdruck gegen eine re 
Wand gefunden: 


\ 


_ 16-0 FEN 
2. ar r ‚far TER 
2 a aofi-b+ ee] 


p) £ -% 

a: 1-22: a DER 7 7 a 5 3 
Ya | 

bzw. 


BL +21 9-31] we” 


Ya 


Es ist zu beachten, daß in Gl. (20) die Größe ö die Gasdichte unmittelbar an der bewegten 
Wand und & die wahrscheinlichste Molekulargeschwindigkeit, wie sie dem Maxwellschen 
Verteilungsgesetz entspricht, an dieser Stelle bezeichnet. Beide Größen werden im allgemeinen 
von den entsprechenden, für das Gasvolumen hinter der Wand gültigen Mittelwerten abweichen 
und ihrerseits noch Funktionen der Wandgeschwindigkeit (v) sein. 


> | Für wird [ 3(0)= 1 | und der Gasdruck auf die ruhende Wand nach GI. (20) 
x 


Bi | 18 a R-T 
EB TE EV er m. ae Due ee we A (22), 


Br: worin R die allgemeine Gaskonstante und 7 die absolute Temperatur ist. 


= f . | Für wird ‚da®o)=1, In . ve Ed = ist. Der Gasdruck p geht. 


also mit waschender Wandgeschwindigkeit asymtotisch gegen Null (p = 0). 
Für d.h, wird 


BN=[3— 2 +) A -DNI- —— —Ro—eNn]... (23) 


Yr:e 


bzw. 


BA) = 3 — 292-1 — O1) —— = 0,05679 —0,31460 : #8. 


Ye 


Für das ideale Gas wird insbesondere | B(l)is = 0,05679 |. Wenn also die Wändge 


schwindigkeit (v) den Wert (&) der wahrscheinlichsten Molekulargeschwindigkeit erreicht, so 
ist der Gasdruck auf die Wand nur noch etwa 6%, des Gasdruckes auf die ruhende Wand unter 


sonst gleichen Bedingungen, insbesondere also bei gleicher Dichte und Temperatur der Gas- 
schicht an der bewegten und an der ruhenden Wand. 


“- “- “- - “- . - “- “ “ - “ “ > 
wm 
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”- 1. Keyartune ER € neck Funktion B Ay 


Es Umformungen und Festsetzungen vor. Wir beschränken uns dabei auf die ein- 


“Um nun die Funktion B(v,) allgemein rechnerisch. auszuwerten, nehmen wir noch einige nah 


‚ 'fachere Form der Funktion nach Gl. (21). Zunächst setzen wir zur APP RUNg en 1-® D(v,) £ 


= Ad; und zerlegen dann den Ausdruck (21) in zwei Teile: 


Bl Et el ae Fu | 


x — Bu( (v, at; 
' Darin ist ’ n 2 
pr Bach 1 x { 
Bu) = ( U+2e): 10 N SE 
7 
‚der Wert von B (v,) für das ideale Gas (0 = 0), während 
DW TR LEBER AR EI N Bade.) 


dem Faktor 92 die Korrekturgröße für die © Abweichung des Gases vom adenjen Zustand 
arstellt. 


Fürik=0.d.1.® = 1 nimmt B(v,) den Wert an: 
2 
Rn . ip 


2 E | 
rn EN EN ET UHREN HE 7.) Bar 
7 
Für zwischen 0 und 1 liegende 92-Werte (0 <92<1) können die zugehörigen B-Werte, 


wie Gl. (25) zeigt, durch einfache Interpolation zwischen den B,-Werten und B,-Werten ae v 
streng berechnet werden. 


Nach der Formel (25) ist die Täbelle 1 der B (u,)-Werte berechnet worden. en Kurven 


von Bild 1 stellen graphisch den Verlauf der Funktion dar. Für einen gegebenen d-Wert wird 
stets für einen bestimmten endlichen o,,-Wert die Funktion B(r,)= 0, während sich für größere 


,  Tabellel. 


0 0 1 1 1 il we 1 


0 1 

0,01 |0,88754 |0,017751 | 0,11172 | 0,79357 | 0,79339 0,79286 0,79197 0,79073 0,77582 
0,04 |0,77730 |0,062184 | 0,21683 | 0,62265 | 0,62203 0,62016 0,61706 0,61271 0,56047 
0,09 |0,67137 |0,120847 | 0,30937 | 0,48235 | 0,48164 |  0,47802 0,47197 | ' 0,46353 0,36200 
0,16 |0,57161 |0,182915 | 0,38462 | 0,36990 | 0,36807 0,36258 0,35344 0,34063 0,18699 
0,25 |0,47950 |0,23975 | 0,43940 | 0,27985 | 0.27745 0,27026 0,25828 0,24150 0,04010 
0,36 |0,39614 |0,28522 | 0,47235 | 0,20901 | 0.20616 0,19760 0,18334 0,16337 | —0,07621 
0,49 |0,32220 |0,31575 | 0,48390 | 0.15406 | 0.15090 0,14143 | : 0,12564 0,10356 | —0,16170 
0,64 |0,25790.. |0,33011 .| 0,47598 | 0,11203 | 0,10873 0,09883 0,08232 0,05921 | —0,21808 
0,81 |0,20309 |0,32900 | 0,45176 | 0,08034 | 0,07705 0,06718 0,05073 0,02770 | —0,24867 
1,0 10,15730 [0,31460 | 0,41511 |'0,05679 | 0,05364 0,04421 0,02848 0,00646 | —0,25781 
1,21 10,11979 |0,28989 | 0,37012 | 0,03956 | 0,03666 0,02796 0,01346 | (—0,00685) | —0,25033 . 
1,44 |0,08969 |[0,25831 | 0,32081 | 0,02719 | 0,02461 0,01686 0,00395 | (—0,01414) | —0,23112 
1,69 |0,06690 10,22308 | 0,27068 | 0,01840 | 0,01677 0,00948 | (—-0,00168) | (—0,01729) | —0,20378 
1,96 |0,04772 |[0,18706 | 0,22252 | 0,01226 | 0,01039 0,00478 | (—-0,00458) | (—-0,01767) | —0,17480 
2,25 |0,03390 |0,15255 | 0,17340 | 0,00805 | 0,00652 0,00195 | (—-0,00568) | (—0,01636) | —0,14450 
4,00 |0,00468 |0,03744 | 0,041333| 0,00079 | 0,00042 | (—0,00071) | (—0,00250) | (—0,00520) | —-0,03665 
6,25 |0,0004070]0,0050875| 0,005446| 0,000049| 0,02 | (—0,00015) BR a —0,005039 


v,-Werte negative _B-Werte ergeben, die keine praktische Bedeutung mehr haben. Es ist von 
Interesse, diese v,,-Werte in Abhängigkeit von Ö zu ermitteln. Nach Gl. (25) kann 'man leicht 
die 9,-Werte berechnen, für die bei gegebenem »,,-Wert die Funktion B{v,) verschwindet, d.h. 
der Gasdruck auf die bewegte a Null wird: 


FERNER: (29). 


2 v2 
% Une rt —Ölo)] 


» 


Die im Nenner. des zweiten Bruches auftretende Größe wre'(i — ®(w,)) hat eine grund- 
legende Bedeutung in der an anderer Stelle 5) entwickelten gaskinelischen Theorie des Aus- 


5) Noch nicht veröffentlicht. 


296 Wand £ 


en 


strömens von Gasen durch Blenden und ist mit X(v,) bezeichnet worden: 


OEL LI EL I En: 


Mit diesem Wert läßt sich Gl. (29) auch in folgender Form schreiben: 


1 1 
2% Yan Xu) 


»=1+ 


000 0% 007 0 WM 2 13 4 SW "VRDUSCUVMBUURDBWUUT M Wu, 
Bild 1 Bild 2 


Nach vorstehender Gleichung ist die Tabelle 2 und die Kurve des Bildes 2für die X (v,)- 
und d,-Werte in Abhängigkeit von ,, berechnet worden. Man kann also mit Hilfe dieser Tabelle 
bzw. Kurve zu einem gegebenen Wert von ® bzw. k sofort denjenigen v,,-Wert finden, für den 
der Gasdruck auf die mit der Geschwindigkeit v, = v,, ' & bewegte Wand zu Null wird, 


Tabelle 2 


lg Ur 
- 1—@(w)] 


0} —o 
0,05 | 0,0010857362 |0,67377089— 2 0,92343156—2 
0,1 0,0043429448 | 0,94818650— 2 0,30110438—1 
0,2 0,0173717798 | 0,19161727— 1 0,45756399—1 
0,3 0,0390865034 |0,30408535— 1 | 0,34317185—1 | 0,59174679—1 
0,4 0,0694871171 |0,35911635— 1 0,42860347—1 | 0,67717841—1 
0,5 0,1085736205. | 0,37975873— 1 | 0,48833235—1 | 0,73690729—1 
0,6 0,1563460135 | 0,37600425— 1| 0,53235026—1 | 0,78092520—1 
0,7 0,2128042961 | 0,35328237— 1 | 0,56608667—1 | 0,81466161—1 
0,8 0,2779484684 |0,31438803— 1 | 0,59233650—1 | 0,84091144—1 
0,9 0,3517785303. | 0,26193490— 1 || 0,61371343—1 | 0,86228837—1 
1,0 0,4342944819 |0,19672654— 1 | 0,63102102—1 | 0,87959596—1 
11 0,5254963231 |0,11983143— 1|| 0,64532775—1 | 0,89390269—1 
12 0,6253840539 | 0,03190594— 1 | 0,65728999—1 | 0,90586493—1 
1,3 0,7339576744 | 0,93343476— 2 | 0,66739243—1 | 0,91596737—1 
1,4 0,8512178453 | 0,82478209— 2 | 0,67599927—1 | 0,92457421—1 
1,5 0,9771625843 |0,70622235— 2 | 0,68338493—1 | 0,93195987—1 
1,6 1,1117938767 | 0,57798065— 2 | 0,68977452—1 | 0,93834946—1 | 1,1525 
1.7 1,2551110527 | 0,44022120— 2 | 0,69533225—1 | 0,94390719—1 | 1,1379 
1,8 1,4071141214 |0,29307736— 2 | 0,70019148—1 | 0,94876642—1 | 1,1252 
1,9 1,5678030797 |0,13666358— 2 | 0,70446666—1 | 0,95304160—1 | 1,1142 
2,0 1,7371779276 |.0,97105678— 3 | 0,70823470—1 | 0,95680964—1 | 1,10406 
2,5 2,7143405119 |0,0075344 — 3 | 0,7218749 —1 | 0,9704498 —1 | 1,07041 . 
3,0 3,9086503371: | 0,82151 — 5| 0,73016 —1 | 0,97873 1 | 1,05020 
3,5 5,3201074033 | 0,4151105 — 6 | 0,7352179 —1 | 0,9837928 —1 | 1,03802 
4,0 6,9487117105 |0,79014& — 8| 0,73885  —1 | 0,987422 —1 | 1,02940 
4,5 8,7944632585 | 0,94702 —10 | 0,74148  —1 | 0,990065 —1 1,02318 


1,024691 


0,75142506—1 


Es sei nochmals besonders darauf aufmerksam gemacht, daß in der vorstehenden Theorie 
der Wert des Korrekturfaktors k 


a. EIITVERRALN a RR 


vr 


für ‚die Abweichung de Bales vom kdealen 7 Pad als nenne von % re ist. 


Dieses trifft in Wirklichkeit nicht zu und zwar selbst dann nicht, wenn die Dichte an der ‚be- 
 wegten Wand konstant bleibt (Z = const), da der Wert von r stark temperaturabhängig istund 


die Temperatur an der bewegten Wand sich im allgemeinen mit v, verändern wird. Wir müssen 


also den k-Wert in unseren Formeln jeweilig durch einen geeigneten Mittelwert für das in Be- 


tracht kommende Intervall der v,-Werte annähern. 


"Die Änderung des k-Wertes mit der Ba und Dichte des Gases folgt den folgenden " Be 


einfachen Gesetzmäßigkeiten:  _ 


Mit zunehmender Temperatur le INS mithin nimmt & mit zu- 
nehmender Temperatur stetls“ab, 


% Bei der Boyle- Temperatur ist insbesondere [a-*=e0] und daher k=1, 9=0|. 


Das Gas folgt also bei dieser Temperatur wiederum dem Gesetz des idealen Gases. 


Unterhalb der Bo yle Temperatur ist: |Jar>eo und k>1; 9>0|. 


Mit zunehmender Dichte (d. i. hneneenden L).nimmt k zu. 


Oberhalb der Bei Temperatur ist [er<eo und k<1; ZZ 


Mit zunehmender Dichte nimmt jetzt kab. 


Der %-Wert ändert sich also mit der Gasdichte stets’so, daß mit zuneh- 
mender Dichte die Abweichungen vom idealen Gaszustand‘ en un 
und mit abnehmender Dichte abnehmen. oe 


Der k-Wert nähert sich dem Wert Null (k= 0) um so mehr an, je Beh die: Temperatur 
und je A E die Gasdichte ist. 


$5. Die Neigung der B Be nel 
Um die Größe des Druckabfalles in Abhängigkeit von der Geschwindigkeit. 4 Be: Wand zu 
untersuchen, bilden wir die Ableitung der Funktion B (v,) nach v,. Aus Gl. (21) folgt 


Eee 
EN ED .e Rai — Ya(1— 9%) - De (1 —D(w)) — 9: ] : 
dv, Yr N 5 
Wir bezeichnen diesen Neigungswert mit tg a 
a N er ER a ee) 


Diese Gleichung hat praktische Bedeutung nur für W-Werte, die kleiner als 9, sind, wo d 
durch Gl. (31) gegeben ist. 

‚Aus vorstehender Gleichung ersieht man sofort, daß der anfängliche Abfall des Druckes 
(v%, = 0) unabhängig von Ö? den Wert 


tg, = — Fn = — 2,257; 20 113295,55,0 
7 
hat. Mit wachsendem v, nimmt die Neigung der B(v,)-Kurve dann ständig ab. Für das ideale 
Gas ist insbesondere: 


en N: ED I en (33), 


m 


während für ein Gas mit der Abweichung 9 vom idealen Zustand eine Korrekturgröße hinzu- 


tritt, die wiederum proportional zu 2 ist. 
E 2 
tgwa=1g8u — - 92 - [Yr - X (u) — v2] TEEN OL 3 (34). 


Schreiben wir den Ausdruck für tga, nach Gl. (33) in der Form: 


re 
er 


are 
a 0,7 Ei a = 2 DO u Ze 


EIk Kr 


so Jäßt er sich mit Hilfe deri in Tabelle 1 Be und sc 
_ Korrekturgröße für die Abweichung ® vom idealen Gaszustand It si 
folgende Form bringen: 


= BE ö . p N: N 
ee ee Er wen: 


= — 2%, [Bo(o) HL — 20) (1 Dw))]- 9 


Mit wachsender Abweichung 9? nimmt die Tangente (tg«) des , al 


nehmende negative Werte an. Die Druckabnahme wird für gleiche ie um so ler“ 


größer die Abweichung (92) vom idealen Gaszustand ist. \ 
Die Tabelle 3 gibt einige nach Gl. (35) berechnete Werte für tar 


 Tabelle3 
Ur 0 0,1 1 6 Bee 0,3 


 tgas | —2,2568 | —1,8794 — 1,5465 — 1,2568 —1,0085 —0,7986 —0,20122 
% 113,9° 118,0° 122,9 | 128,5° 143,7° 141,4° 168,6° 


$ 6. Der Gasdruck auf die bewegte Wand bei Kompression des Gases 

Die Gleichungen (20) und (21) liefern auch sofort die Gesetzmäßigkeit für die Zunahme 
des Gasdrucks auf die bewegte Wänd, wenn diese das Gasvolumen komprimiert. Man braucht 
in diesen Gleichungen nur v, durch seinen negativen Wert (—v,) zu ersetzen; denn die Integrale 
der Gl. (14) sind in diesem Falle von —v, bis +00 zu erstrecken, da auch noch die Moleküle 
gegen die Wand stoßen, welche eine von der Wand fortgerichtete Geschwindigkeitskomponente (z) 
haben, deren absolute Beträge kleiner als derjenige der Wandgeschwindigkeit v ist [|z] sv]. 
Man erhält daher jetzt für den Gasdruck die Gleichung: 


. (37), 


worin wiederum v%, = v:& und v der absolute Betrag der Wandgeschwindigkeit ist. 
Für den idealen Gaszustand vereinfacht sich die 2. Bo zu der Form: 


Be =L FO +2 + N LEE A 


Für a) ergibt sich daraus wiederum der Wert B(0) = 1. 


Tabelle 4 
R » 
Yr ++] 1+9(w) | DEN r 2v (1 + Dr) 2(1+2v% 
rt 
0 1 1 

0,1 1,02 1,11246 0,11172 0,0222492 1,24643 2,04 
0,2 1,08 1,22270 0,21683 0,097816 1,53735 2,16 
0,3 1,18 1,32863 0,30937 0,239153 1,87715 2,36 
0,4 1,32 1,42830 0,38462 0,457085 2,27009 2,64 
0,5 1,50 1,52050 0,43940 0,76025 2,72015 3,00 
0,6 1,72 1,60386 0,47235 1,15478 3,23099 3,44 
0,7 1,98 1,67780 0,48390 1,64424 3,80594 3,96 
0,8 2,28 1,74210 0,47598 2,22989 4,44797 4,56 
0,9 2,62 1,79691 0,45176 2,91099 5,15966 5,24 
1,0 3,00: 1,84270 0,41511 3,68540 5,94321 6,00 
1,1 3,42 1,88021 0,37012 4,55011 6,80044 6,84 
1,2 3,88 1,91031 0,32081 5,50169 7,73281 7,76 
1,3 4,38 1,93401 0,27068 . 6,53695 8,74164 8,76 
1,4 4,92 1,95228, 0,22252 7,65294 9,82774 9,84 
1,5 5,50 1,96610 0,17840 8,84745 10,99195 11 
2,0 9,00 1,9953223 0,041333 | 15,962579 17,999234 18 
25 | 135 1,9995930 - 0,005446 24,9949125 26,999952 27 


a En 3 at Sn ud ul Zub cn a 0 a 


Für o Re) geht jetzt B, ) mit 1 +2. v Bogen unendlich, Für Binänglich role tr & 


Werte (2 1,5) wird jetzt allgemein mit großer Annäherung 


h B-W-2U +20 WEld. 2.0: 3. 


Für den speziellen Wert ?, = 1,5 wird der Nr der vereinfachten Berechnung nach vor- 
stehender Formel bereits kleiner als 19/50, der genaue Wert beträgt B(— 1,5) = 10,99195, während 
sich aus der Gl. (39) berechnet B(—1,5)= 11. Für «= 1 bleibt der Fehler noch unter 1%. 


Voraussetzung für die Berechnung des Druckes nach der Gl. (37) ist wiederum, daß de 


Größen ö und a2, d. h. die Dichte und Temperatur des Gases unmittelbar an der bewegten Wand 
in jedem Augenblick gegeben sind. Allgemein kann man die von unserer "Theorie gelieferte 


Funktion B(v,) in ihrer Bedeutung folgendermaßen klarstellen: 


Wenn eine mit der Geschwindigkeit v be- gu, 
wegte Wand, durch die das Gas entweder ex- ®7 

pandiert oder komprimiert wird, plötzlich fest- 
gehalten wird, ihre Geschwindigkeit also mo- 
mentan zu Null wird, so ändert sich der Gas- 
druck auf die Wand im gleichen Moment im Ver- 
hältnis B(e,) :1. Im nächsten Augenblick wird 
sich dann dieser Druck natürlich im allgemeinen 
noch weiterhin ändern, da infolge der im Gas- 
raum vorhandenen : Strömungsvorgänge sich die 
Dichte (ö) und die Temperatur (a2) an der Wand 
noch verändern werden, und zwar erfolgt diese 
Änderung nach einer grundsätzlich anderen Ge- 
setzmäßigkeit, wie sie durch die Funktion B (v,) 
dargestellt wird. 


Tabelle 4 enthält die nach Gl. (38) berechneten Werte für. B,(— ©,). Die für die Abweichung 
vom idealen Gaszustand erforderliche Korrektur der B,-Werte beträgt: 


AB= Bi) — Bi— 9) = — 2ur [1 + D(u)] 92 


Bild 3 


Für den Wert #=1] geht B(— v,) schnell gegen den Grenzwert 2, nachdem es für v, = 1,0 
das Maximum B, (— 1) = 2,25781 erreicht hat. Bild 3 zeigt die zugehörigen Kurven. 


$7. Zum Widerstand bewegter Flächen 

Ich weise noch auf die besondere Bedeutung der vorstehenden Theorie für den Widerstand 
der in Gasen bewegten Flächen insbesondere auch für fliegende Geschosse hin, denn bei letzteren 
treten Werte von ®, auf, die beträchtlich über die Einheit hinausgehen können. Die Theorie 
zeigt, daß das gebräuchliche quadratische Luftwiderstandsgesetz offenbar nur ein idealisiertes 
Grenzgesetz ist, und daß mit besserer Annäherung für sehr hohe Geschwindigkeiten der Wider- 
stand proportional zul +2(1 — 0°) » v7 gesetzt werden kann. 

Allgemein erhält man den Widerstand W(v,) einer in einem Gase bewegten Fläche (Kör- 
pers) durch die Differenz der Gasdrucke auf die Vorder- und Hinterseite der Fläche. Bezeichnen 
wir in den obigen Gleichungen die Größen, die sich auf das komprimierte Gas an der Vorderseite 
beziehen, mit dem Index c und diejenigen Größen, die sich auf das expandierte Gas an der Rück- 
seite der Fläche beziehen, mit dem Index e, so läßt sich der Widerstand W(v,) der bewegten 
Fläche in Abhängigkeit von seiner Geschwindigkeit (v,) allgemein darstellen durch die Differenz 
der beiden Gleichungen (37) und (20). Wir setzen darin noch: 


We)=t 9 = Me’ Bi-U)— me’ Bü) -.. .. (42), 
worin 


2 
B—t)= [1 +2(1- 9) u] (1 +0(w)) er 
3 


Bw) =[11+2(1-—- 9) :v]- Us) ne x 


ist. 


SE 2 UBS ETUW UN ZRUKLUSULNT 1 WU . 


,.. "und daher: 


 Ausdiesen Gleichungen kann man Er En zgesetze für 
zwar einmal bei sehr kleinen Geschwindigkeiten (w<1 wenn die 
a sten Wand sowohl auf der Vorder- v8 auch auf der R F 


Es ist aber: 


B(— vw) — By) = og +2(1— 8) v7] - a = rn og 


7 Ku 


Die zweite an ergibt sich für hinlänglich große ER (% > >> re 
_ da dann der Druck auf die Rückseite der Fläche De dem Druck auf die Vorderseite ver- - 


nachlässigt und gleich Null gesetzt werden kann. . Dann nimmt die Widerstands- 


funktion folgende Form an. j 


Wu) =. 9.B(— v) = Tre 1 E ©(r,)) | +21 —H)v)+ ed wi 


Yr 


Da darin © (v,) Ar nahe gleich der Einheit geworden ist und v, SE ? bereits kleine Werte 


Hat, so kann praktisch noch einfacher geschrieben werden 


ET A BEE ee (47). 


. Dabei ist zu beachten, daß der Wert von 92 infolge der starken Kompression der Gas- 
schicht vor der bewegten Fläche bei sehr hohen (Überschall-)Geschwindigkeiten größere Werte, 
die sich der Einheit annähern, erreichen kann, so daß der Widerstand wesentlich langsamer, als 
es dem Quadrate der Geschwindigkeit (v}) entspricht, zunehmen muß. 

Für das zwischenliegende Intervall der v,-Werte (0,01 <v, < 1) ist ein genaueres Wider- 
standsgesetz nur unter Berücksichtigung der Änderung der Gasdichte (6) und der Temperatur 
in der Gasschicht an der bewegten Fläche zu gewinnen. 


Anhang 
Die vorstehend entwickelte Theorie liefert nebenbei noch eine sehr einfache und zweck- 
mäßige Näherungsformel zur Berechnung des Gaußschen Fehlerintegrals ®(z) für große 
Werte des Arguments. Wie Tabelle 1 zeigt, nimmt der Wert der Funktion B (v,) mit wachsendem 
Argument sehr schnell ab und ist bereits für v, = 2,5 kleiner als 10. Daher erhält man aus 
Gl. (26) eine sehr einfache Gleichung zur Berechnung ter ®-Funktion für größere Argumente, 


wenn man darin setzt. So folgt‘ 


Nach dieser Gleichung sind einige P(x x)-Werte berechnet und in der nachstehenden Ta- 


belle 5 zum Vergleich mit Werten, die einer Funktionstabelle der ®-Werte entnommen sind, 
zusammengestellt. 


Tabelle5 


1-9 

a 4,593 - 10° | 4,034 10”4 | 2,199 - 10-° | 7,411 - 10-” | 1,539 - 10-8 | 1,964 - 10-10 | 1,5364 - 10-% 
1—® 

nach Tabelle) 4,678 10° | 4,070 10-4 | 2,21-10-° | 7,431. 10-7 | 1,542 - 10-® u 


Eingegangen: 24.11. 1948, 


bisherigen theoretischen Arbeiten über Propellerschall — auf die mathematische Behandlung 
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Der Propeller als rotierende Schallquelle 
Von H. MerbtundH. Billingin Göttingen 


„Es werden M elhoden entwickelt zur Berechnung des Schallfeldes rotierender Luftschrauben in Abhängig- 
keit von den geometrischen Abmessungen, Drehzahl, Anstellwinkel und Forischritisgrad der Luftschrauben. 
Die theoretischen Ergebnisse werden mit Ernsthausens Messungen an Modellufischrauben verglichen. 

Methods. are developed to caleulate the field of sound of rotaling air screws dependent on the geometri, 


data, number of revolution, angle of incidence and ratio of advance of the screws. T’he theorelical results i 


„are compared wilh Ernsthausen’s measurements on model screws.. _ 

On d£veloppe des meihodes pour caleuler le champ de son d’helices en je des dates geom£lriques, 
du nombre de r&volutions, de l’angle d’incidence ei de l’argument de similitude des helices. Les r£sultats 
ihoreliques sont compares avec les mesures X’Ernsthausen faites ü des mod2les d’helices. 

‚PaspaÖatssamrca MeToxsI pacyeTa 3BYKOBOTO MONA BPAINAMINNUXCH BOSAYINHEIX BHHTOB 
B 34BHCHMOCTU OT TEOMeTPHYeCKHUX MAHHBIX, OT YUCAHA O00POTOB, OT YTAIA ATaku u ara 


suHTa, Teoperuyeckue pesyIBTaTsl CPABHHBAWTCH C PesyıbTaTaMH ucukTaunl OPHeT- 


Tays3eHa Ha MONeIAX BO3AYIHHEIX- BEHTOB. 


1. Problemstellung und Differentialgleichung 


| Die folgenden Ausführungen sind im wesentlichen in den letzten Kriegsjahren erarbeitet 
und bisher lediglich in einer Reihe von noch unzugänglichen Berichten niedergelegt [1]. 


Zahlreiche Messungen an Flugzeugen und Modellschrauben haben gezeigt, daß sich der 
Propellerschall im wesentlichen aus zwei separaten Anteilen zusammensetzt, dem durch Wirbel- 
ablösung am Schraubenblatt erzeugten Drehgeräusch und dem sogenannten Drehklang. Dieser 
besteht aus einer Reihe von harmonischen Teiltönen und hat seine Ursache in einem mit dem 
Propeller umlaufenden Druckfeld. Da das Wirbelgeräusch wegen seiner geringen Stärke und 
seiner hohen Frequenzen in der Luft schnell absorbiert wird, beschränken wir uns — wie die 


des. Drehklanges. 


Zwei Gründe für die Entstehung des Propellerklanges sieht man unmittelbar, wenn man. 


ein Raumelement betrachtet, durch welches ein Propellerblatt hindurchschlägt. Einmal wird 
durch das Propellerblatt die Luft verdrängt, und zwar nach beiden Seiten gleichmäßig. Diese 
Verdrängungswirkung hängt praktisch nur von der Dicke und Tiefe des Propellerprofils ab und 
führt akustisch gesprochen zu Nullstrahlern. Zweitens erteilt der Vortriebspropeller dem Luft- 
element einen Impuls in Richtung der Propellerachse. Dieser gerichtete Impuls, der jedoch nur 
bei belasteten Propellern auftritt, führt akustisch zu Strahlern erster Ordnung. Es .erweist sich 
nun. als sehr nützlich, auch in der mathematischen Behandlung die Verdrängungswirkung und 


die Schubwirkung gesondert zu berechnen und dementsprechend das Schallfeld aus der Wirkung - 


eines unangestellten dicken und eines angestellten unendlich dünnen Propellerblattes zusammen- 

zusetzen. Der mathematische Grund zu dieser Aufteilung ist unmittelbar klar: Um den-un- 

angestellten dicken Propellerflügel bildet sich keine Zirkulation aus, wenn man von der das 

Wirbelgeräusch erzeugenden relativ schwachen Verwirbelung absieht. Hier gilt mithin mit Aus- 

nahme des Propellerortes im ganzen Raum für das Geschwindigkeitspotential 9 die homogene 
Wellengleichung. 

1 29 
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Hinter dem angestellten Propeller hingegen ist im gesamten Bereich des Propellerstrahles 

rot vo #0. Es muß daher mit Druckpotentialen gerechnet werden, wie vom zweitgenannten 

Verfasser in [2] ausgeführt wurde. Um Druckpotentiale anwenden zu können, soll die Gültigkeit 

der Wellengleichung (1) für den Schalldruck im gesamten Raum um den angestellten Propeller 

zunächst vorausgesetzt werden 


RE N NE El 
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Gl. (2) ist jedoch gemäß [2] nur dann im ganzen Raum erfüllt, wenn die Felder der vom Propeller 
erzeugten freien Wirbel zeitlich stationär und raumfest angesetzt werden dürfen. In der Praxis 
sind beide Voraussetzungen nicht exakt erfüllt, da die im Propellerstrahl auftretenden Wirbel 
zusätzlich mit der Strahlgeschwindigkeit nach hinten abschwimmen und sich auflösen. Die 
Strahlgeschwindigkeiten sind jedoch im Fall des fliegenden Flugzeuges stets klein im Ver- 
geich mit der Fluggeschwindigkeit, außerdem nimmt die Zusatzgeschwindigkeit im Propeller- 


- strahl mit anwachsender Vorwärtsgeschwindigkeit (bei konstant gehaltenem' Schub) ab. Bei 


der Berechnung des Schallfeldes vorwärts bewegter Propeller kann deshalb von der zusätzlichen 
Bewegung der Wirbel im Strahl abgesehen werden. Die Gültigkeit der Gl. (2) läßt sich dann 
für den ganzen Raum voraussetzen. Be 
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"nachlässigung der Zusatzgeschwin or im Strahl zu erheblichen Fehlern, nase cn... » 
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g, Dar Propslier als rotierende 


Bei der Schallfeldberechnung eines im 'Stand laufenden Propellers führt dagegen die Ver- 


x h 
k 


‚auch in diesem Fall mit Hilfe einer 
wie später gezeigt wird. MS N PR TEER SEN 

Zur Berechnung des Schallfeldes einer Luftschraube wird die Singularitätenmethode ver- 
wendet. Dazu müssen zunächst kreislinig bzw. schraubenlinig umlaufende Singularitäten kon- 
stanter Ergiebigkeit berechnet werden, um mit ihnen die Randbedingungen am Propellerblätt 
zu erfüllen. m 


eorie höherer Näherung brauchbare Ergebnisse 


“ | x Für den unangestellten Propeller wird als Randbedingung die Normalgeschwindigkeit am 


Propellerblatt gewählt, Es genügt dabei, in der folgenden Rechnung das genaue Flügelprofil 
durch einen einfachen Quell-Senkenkörper von gleicher Dicke und Tiefe zu ersetzen. Eine gegen- 


a seitige Beeinflussung der unangestellten Propellerblätter bleibt bei dieser Näherungsrechnung 


unberücksichtigt. | fi a RE 
Als Randbedingung für den angestellten Propeller ebenfalls die Normalgeschwindigkeit zu 
nehmen und diese Randbedingung wie in [2] mittels Drucksingularitäten zu befriedigen, ist bisher 
noch nicht versucht worden und erscheint recht schwierig. Leichter gangbar erweist sich der 
folgende Weg: an 
Mittels der jetzt gut ausgearbeiteten Propellertheorie errechnet man sich unter Benutzung 
der Wirbeltheorie die Wirbelbelegung und damit die Auftriebsverteilung längs des gesamten Pro- 
pellerblattes. Damit steht als Randbedingung der Drucksprung am Propellerblatt zur Verfügung. 
Zur Berechnung des Schallfeldes braucht man dann lediglich die nach der Propellertheorie be- 
rechnete Wirbelbelegung durch eine proportionale Belegung mit Drucksingularitäten zu ersetzen 
und über die Belegung zu integrieren (vgl. [2]). Die Drucksingularitäten laufen dabei mit dem 


‘ Propellerblatt um, Die Genauigkeit der Schallfeldberechnung hängt für den angestellten vor- 


. wärtsbewegten Propeller dann im wesentlichen nur von der Güte der Rechenmethode ab, mit 


deren Hilfe die Auftriebsverteilung am Propellerblatt bestimmt wird. Für den akustischen Teil 
der Rechnung wird als einzige zusätzliche Voraussetzung die Gültigkeit der Wellengleichung (2) 
im gesamten Raum: außerhalb der Propellerblätter benötigt. 

Gl. (1) und (2) werden sicherlich ungültig beim Auftreten von Verdichtungsstößen. Eine 
etwaige Anwendung der im folgenden abgeleiteten‘ Formeln auf Propeller, deren Spitzen mit 
‚Überschallgeschwindigkeiten umlaufen, erscheint deshalb gefährlich und ist noch nicht näher 
diskutiert. 


2. Kreis- und schraubenlinig bewegte Singularitäten 

Zur Berechnung des Feldes kreis- bzw. schraubenlinig bewegter Punktsingularitäten sind 
verschiedene Ansätze versucht worden. Eine Ausnützung der Invarianz der Wellengleichung 
gegenüber der Transformationsgruppe der konformen Abbildung des vierdimensionalen Raumes, 
welcher Weg sich in [2] für die geradlinig gleichförmige Bewegung als so besonders elegant er- 
wiesen hat, ist nicht möglich. Es steht für die Kreisbewegung keine geeignete Transformation 
zur Verfügung. 

Ein für die Kreisbewegung zum Ziel führender Ansatz von Ernsthausen [3] belegt 
die Kreisbahn dicht mit ruhenden Elementarstrahlern, die jeweils nacheinander beim Durchgang 
der bewegten Singularität einen kurzen endlichen Störimpuls aussenden. Die Größe und Dauer 
des Impulses richtet sich nach den vorgegebenen Randbedingungen, die erfüllt werden sollen. 
Man erhält eine Aufspaltung in harmonische Teilschwingungen, wenn man den zeitlichen Stö- 
rungsverlauf am Ort eines Elementarstrahlers in Form eines Fourier-Spektrums darstellt. Für 
die Berechnung des Potentials in einem beliebigen Aufpunkt müssen dann die Wirkungen aller 
Einzelstrahler phasenrichtig addiert werden. 

Wir wollen diesen anschaulichen, aber mathematisch etwas umständlichen Weg hier nicht 


weiter verfolgen, sondern für eine zunächst ganz allgemeine Quellbewegung ausgehen von der 
klassischen Wellengleichung d’Alemberts 


19%y ’ 
A V— = 9: =g9 (z, Y,2, t) ER a a a eh ae ee. (3). 


Darin ist x ein beliebiges Potential und g (x, %, 2, t) bedeutet die.die Störung verursachende Quell- 
verteilung. Eine Lösung von (3) ist nach Poisson 


de 
ya | —- dl... un (4). 
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Zeit nn zn integrieren. a über Ken ann zur retardierten Zeit von Quellen besetzten 


' Raum. Die Integration ist im allgemeinen sehr schwierig, da Lösungen gesucht werden, die den € 
Quellort: und das gesamte Feld der Störung zur Aufpunktszeit t angeben. , 

' Für eine Quelle konstanter Ergiebigkeit, welche auf einer beliebigen, in sich BeschienE < 

- Bahn mit gleichförmiger Geschwindigkeit umläuft, hat Schott [4] eine elegante Umformung 


der Poissonschen Lösung (4) angegeben, welche von Küssner [5] für die Theorie ds 
Propellerschalles. nutzbar gemacht Va ist. Schott entwickelt in der Poisson schen BR; 
Lösung für die mit. der Periode Da ee bewegte Singularität die Quellverteilung E 
ry e 
i o( ET in eine Fourierreihe zu. > 
R | o( 2) = PH a ar NE nei (5) 
A, > Mm=— 00 fi 
Dabei ist . hr 
RA Ä 
le wi 
oder mit T=t— = unter Ausnutzung der Periodizität vn gay, . 
e f { q i ? h Be ee: 
a — imo 7+ he | Ä, u, a ER Be: 
Sr | I PR en (7). a - 
' Damit wird (5) | | ; En 
r' T 1 = imo (1 2) ji 
o( ee e NE) AT IT re N (8) 
Die Integration über die retardierten Zustände ist dabei ersetzt durch eine. Integration Be 
über eine ganz Periode der Bewegung. Beachtet man weiter, daß die gesamte Quellstärke kon- > 


stant bleibt, so folgt schließlich aus (4) und (8) das Potential für eine periodisch mit der Periode 7 
bewegte Quelle konstanter Quellergiebigkeit C 


ed | Ne 
nd es in AUREAC REN NE (9). E. 


mM= > % 


Der Strich bei &’ bedeutet, daß für m — 0) nur die Hälfte des Wertes einzusetzen ist; m gibt die a 
Ordnungszahl der Schwingung; r’ ist der Abstand zwischen Aufpunkt und Quellpunkt zur Zeit r ö 
und für die Kreisbewegung leicht in Abhängigkeit von 7 anzugeben. Dazu sind mit Küs S- Be 
ner[5] der Aufpunkt und der  Quellpunkt in Zylinderkoordinaten gegeben zu ER 
2". y=ssiny .2=8c0sp \ (10) ri 
x y=Rsin(wor +6) 2= Rcos(wr-+öÖ) I ER : E 

"Das Quadrat des Abstandes zwischen diesen beiden Punkten ist 


r2=r% —2Rscos (wor +d5.— y) | 


EEE N SS N RE CHEND NER MN, NE (11) 
3=(2— a ++ R? | 
Führt man noch die dimensionslosen Parameter 
EEE AR BE ee (12) 
C rs 


ein und setzt or +6 — y =, so folgt aus (9) als retardiertes Potential für eine mit der Win- 
kelgeschwindigkeit auf einem Kreis vom Radius R umlaufende Punktquelle 


„= >" hm(a,0)explim(wE H+ö—Yy] -» -» re. (13). 


52 re 


304 Merbt u. Billing, Der Propeller als rotierende Schallquelle 


Dabei ist 


exp [im (—L—afl- ne (14). 


f 2 Ä d 1 
T TTUM MR T 
nr) - | 
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Für die Propellertheorie benötigt man das Potential für z äquidistant auf derselben Kreis- 
bahn bewegte Quellen. Man setzt dann in (13) für den Phasenwinkel der k-ten Quelle 


= ty helm... are 


und summiert über k. Dabei ist, wie man leicht einsieht, 
| 2 Over 0° für m +0 (mod) 
2, pP aa e| je | f = i 
k=1 a 2 ür m= 0 (mod z) 


Damit folgt aus (13) für das Potential von z Quellen, die sich äquidistant auf derselben Kreis- 
bahn bewegen ; 


Kr e be Amz(o,0)expfemz(wt Hö—y)- -...-.- (17). 


z IST rg Ken 


Für die Funktionen Am; sind in [6] Reihenentwicklungen angegeben worden, die eine Berech- 
nung der Potentiale gestatten. : 

Betrachtet man weit entfernte Aufpunkte (Abstand > 5 Bahndurchmesser), so können 
wir in (14) näherungsweise setzen 


im Zähler yl—2ocosö =1—ocosZ 
im Nenner Yl—2ocosc—1 


8 s 2—& z—-2 
== — = 0509, re 


ER = sin® 
15 Li Tg 


und, da Re<(s—-z” +32=r ist, folgt 
(© = Breitenwinkel des Aufpunktes). 


Führt man noch die Mach sche Zahl M = = ein, so erhält man als Näherungslösung 
für das Fernfeld von 2 kreislinig,bewegten Punktquellen 


08: 


KT ar 


ID ir: Im: (mzMp cos ©) exp imzo ((— 2 ei, + (18) 


- MU na 
(Jmz = Besselsche Funktionen). j 

Für die Propellertheorie benötigt man weiter Lösungen für kreislinig bewegte Quell- 
Senken-Anordnungen und Dipole. Man erhält sie in einfacher Weise durch Differentiationen 
aus den bereits abgeleiteten Formeln (17) und (18). Ordnet man den z umlaufenden Quellen z 
um den Phasenwinkel d ö verschobene Senken zu und geht zum Dipol über, so wird das Fernfeld 
für z solche auf einer Kreisbahn umlaufende Quell-Senken-Anordnungen 
20,dd < na n pi ( r ] 

Dmi Jmz (mz Mpx cos ©) exp |limzw(|t e En . (19). 


m—=1 


_ 9% ee. 
Ma DER Dr 


Ebenso ergibt sich für z kreislinig umlaufende Dipole mit ihrer Achse senkrecht zur Anströmung 
in Richtung der x-Achse das Potential 


0) 

Ka = an da, = 

Aa0,deo. sa r  ö6+ 

ern sin © 2 Ma? ms (mz MR cos 9) exp imeo ( eh a] (20) 
(quadratisch mit dem Abstand abnehmende Glieder sind hierin vernachlässigt). 

Um schließlich aus der Kreisbewegung Lösungen für schraubenlinig bewegte Singularitäten 
zu erhalten, braucht man lediglich gemäß [2] die Gleichungen (17) und (18) einer Lorentz- 
Transformation zu unterwerfen. Wendet man die dort für den ruhenden Beobachter angegebenen 
Transformationsformeln 


= (21) 


Vx 
’ [2 ’ 7 v f 
= a y- Y 2 F c? (2) 


F t= —— — mit ?=-1— 
% X F: 2 it % 1 


Br 


ae ix Laer 
Br, 


Ri, 


ui = nm — V,t 


& d. 2. € Quelle regte: mit der Winkelgeschwindigkeit @ oe einem Kreis mit deal Radius E 
 B=R:% und gleichzeitig mit der ‚Geschwindigkeit. V, in der Richtung der negativen x-Achse. 


i Der: Fortschrittsgrad der Quelle beträgt = ER 
'Lorentz- Transformation für den mitbewegten Beobachter 


SB: 


RE A wen XM . Le 


a A et mit ut RB HR 


Muh Sy, ER RE 0%? 


ersibteh sich aus (18) das Potential für z Sehraubehlmie bewegte Quellen es Ergiebigkeit C-% 
. bei re a ee Zu 0% 


‚ C'2 


at Imr cos ©’ Ytg? 0’ zur 


> m z mzM H M(«' nr 128, 24-425 r3 de 
> D 2 Felt ken Ba EN er 3% ut (24). 
Sa 


Dabei ist Mx die Machsche Zahl der Projektion der Quelle in die y, z-Ebene, M= — = 


m=0 


: die Machsche Zahl der Quellprojektion auf die z-Achse, #? = 1 — M?. 


Wie man leicht einsieht, gelten die Formeln — ebenso wie alle bisher abgeleiteten - — sowohl 
für umlaufende Singularitäten des Geschwindigkeits- als auch des Beschleunigungspotentials. 
Sie sind also zur Darstellung des nicht angestellten Propellers wie auch. des angestellten. Pro- 


pellers geeignet. 


3. Drehklang von Luftschrauben 
a) Der nicht angestellte Propeller im Standlauf 


Um die Strömung um den nicht angestellten Propeller nachzubilden, müßte man den Pro- 
peller in solcher Weise mit kreislinig umlaufenden Quellen und Senken des Geschwindigkeits- 
potentials besetzen, daß an den Propellerblättern die Randbedingung für. die Normalgeschwin- 


digkeit erfüllt wird. Als ausreichende Näherung erweist es sich jedoch, die Quellsenkenverteilung - 
. durch eine Belegung des Propellerflügelrandes mit einer Quell- und einer gleichstarken Senken- 


linie zu ersetzen. Der Winkelabstand zwischen diesen beiden Linien beträgt beim Radius 
R: & =dö(l (R)= Tiefe des Propellerblattes). Die Intensität C|der Quellbelegung beim Radius R 


ee man aus der Dicke D(R) des Propellerprofils bei Rechnung in ebener Näherung (Strei- 
fentheorie)!) zu 


Schließlich kann man bei nicht zu tiefen Profilen für die ersten Harmonischen unter Beibehaltung 
des Momentes C, dö = Diw die Quell- und Senkenlinie zusammenfallen lassen, also zur Dipol- 
bildung übergehen und erhält durch Einsetzen von Gl. (25) in Gl. (19) und Integration über den 
Radius.als Geschwindigkeitspotential für das Feld eines ya nicht angestellten Propellers 


en exp » [ine(t— [m (m2MpcosO)IDdR . (26). 


7 Er mr 
Setzt man D= Dund ds — (d= mittlere Flügeltiefe) über den Propellerradius konstant und 
substituiert in (26) . 5 | n 
2 m: 00:0 = me” 000, 


so vereinfacht sich die Gleichung (26) zu i 

6 — 1 er 
gu Ze Din, I minetiesp[inzol— + Ze je, ren 
’ 9 


& 


1) Vgl. z.B. Prandtl: Strömungslehre. 
20 


0 R = e” RR ER cost =1, d.h. or +0= v a, A EN a 


' Mit‘ der in [2] EBEN | 


RE I N (5). 


wire 


m Standlauf(l.] 


er 


RER; NT A 
b)Der angestellte Propellerim dlauf(1.Näh ! 
In den folgenden Abschnitten wird die Gültigkeit der Wellengleichung (2) ür 
Raum außerhalb der Propellerblätter vorausgesetzt, der Einfluß der Zusatzgeschwindig 
.- im Schraubenstrahl also vernachlässigt. Für den im Stand laufenden Propeller fällt d 
nachlässigung durchaus ins Gewicht. Gleichung 29), (30) und (33) sind deshalb nur als 
äherung zu bezeichnen. eg EI wi Dr 
e Um unter diesen Voraussetzungen das Schallfeld eines angestellten Propellers zu berech- 


. » - were 
nen, ermittelt man mittels der Propellertheorie zunächst die Zirkulationsverteilung über den az 


Propellerradius. Darauf belegt man den Propellerradius so mit Druckdipolen (Dipolachse in 
Richtung der Propellerachse, d.h. x-Achse), daß die Querkraft auf ein Element der Dipollinie 


gleich der bekannten Querkraft auf das Propellerelement L= o VI'(R) ist. Aus dieser Forde- TR 


. rung folgert, wie man leicht überlegt, als Wert für die Dipolkonstante c (R) -- . 
| L(R\AR=C(RBJdÄR=gVTdR. . .. 22.2... (28). 


Das Schallfeld des angestellten Propellers ergibt sich dann unmittelbar, wenn man C{R) aus wo 


-Gl. (28) statt O,da, in Gl. (20) einsetzt und über den Propellerradius integriert. x hat dabei 
die Bedeutung des Schalldruckes p. Man erhält für den Schalldruck in großer Entfernung vom 
Propeller - 


Rg F 
1 = ch sin ® 3 mist | Jmz(mzMpzcos ©) T(R) RdR- 
ER ‘ i 


a 


I. 


| Setzt man in erster Näherung für die Auftriebsverteilung I’ = z loR (%, = mittlerer Auftriebs- 
beiwert) und substituiert wie in (26), so vereinfacht sich die Gleichung (29) zu 


£ 
Per I all 02 mimH 3 | Amtstazıcnp imzw Kart To (30). % 


Die Integrale in Gleichung (26) und (29) müssen im Spezialfall numerisch ausgewertet werden, 
während die Integrale in Gl. (27) und (30). 


D,,m:(d) = f u FFIC)L E Speer en (31) 


von Ernsthausen [7] berechnet und für den interessierenden Bereich tabuliert sind. 

Zur Bestimmung des gesamten Schallfeldes eines angestellten dicken Propellers überlagert 
man in schon besprochener Weise das Feld eines angestellten, aber unendlich dünnen Profils dem 
Feld eines unangestellten dicken Propellerprofilst). Zur Berechnung des Fernfeldes hat man also 
lediglich die Teillösungen (26) bzw. (27) und (29) bzw. (30) zu superponieren. Zur Ermittlung des 


resultierenden Schalldruckes muß man mittels der Beziehung p = o a dazu in Gl. (26) zunächst 


‘ das Geschwindigkeitspotential in den Schalldruck umrechnen. Das Druckfeld des „dicken‘ 
n 
nr 


ben. Für den Schalldruck p= /p; +} des angestellten dicken Propellers im Standlauf erhält - 
man dann für das Fernfeld . 


und das des „‚angestellten‘ Propellers ergeben sich gegeneinander um die: Phase — verscho- 


nten imo N 
p= A exp |imzw denen ur; 


[02] 
Rg 2 n Rs 2 (32). 
sin © 
. \f Didmsmeitncca an) + (mie [ Tim Ran. meitacn ar) 
ö ö 


') Es liegt der Einwand nahe, daß durch die Anstellung des Propellers das durch die Dickenwirkung er- 
zeugte Schallfeld verändert wird, also keine lineare Superposition gestattet ist. Dieser Einwand läßt sich 
damit begründen, daß im Propellerstrahl jetztrot v>£ 0 ist und mithin das von der Pro llerdicke herrührende 
Schallfeld nicht mehr mittels Geschwindigkeitssingularitäten berechnet werden darf. Zur Vermeidung dieser 
Schwierigkeit kann man das durch Dickenwirkung am Propellerblatt erzeugte Druckfeld als Randbedingung 

‚Wählen und mit Druckpotentialen aufbauen. Damit verzichtet man auf Wellenglaickuns (1). Die auf diesem 
"Weg gewonnenen Endformeln unterscheiden sich nicht von den im folgenden abgeleiteten. 


Lei. Die Zn; 


APR FREMDER 


ne 
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ke, Morde Biliag, mern Propeller als rotierende Schallqulie 


er Da die Schallwellen für. das Fernfeld als eben angesehen ach können, folgt aus. Gl. (32) die 
'Schallschnelle o mittels der Beziehung p= oc», womit sich auch die Schalleistung des Pro- 

pellers unmittelbar berechnen läßt. 

Eine erste Theorie des Propellerschalls ist in einer Reihe von Arköten von Ernst- 

hausen [3], [7], [8], [9] auf einem anderen Wege entwickelt worden. Statt der ae ; 

. methode benutzt er die Kirchhoff a Lösung der Wellengleichung 


exp [io ( (2) | exp N, 
d c 1 0 | 
a, dS +; H ————Jäs . (83). 


Die ration wird über eine. die Schallquelle eng einschließende Kontrollfläche ausge- 


führt. 2% und x, sind dabei die Randwerte der Geschwindigkeit und des Druckes am Schrauben- 


kreis. Das erste Integral enstpricht unserem Ansatz für den nicht angestellten Propeller, das 
zweite dem für den angestellten Propeller. Da die notwendige Unterscheidung in der Anwendung 
des Geschwindigkeits- und Beschleunigungspotentials damals unbekannt war, ist in [3] die Be- 
stimmung der Randbedingungen ungenau. Zum Angleich an die Formeln (27) und (30) sind die 
Endformeln Ernsthausens deshalb für den unangestellten Propeller mit dem Faktor im 


und für den angestellten Propeller mit a zu multiplizieren. 


n 


c) Der vorwärts bewegte Propeller 


In den praktischen Anwendungen interessiert vor allem das Schallfeld des Propellers, der 
sich längs seiner Achse vorwärts bewegt. Die Vorwärtsbewegung verursacht bei den im modernen 
Flugbetrieb vorkommenden Geschwindigkeiten bereits eine recht beträchtliche Änderung des 
Propellerschallfeldes. Die Anwendung der für den Standlauf abgeleiteten Formeln auf das Schall- 
feld eines vorwärts bewegten Propellers ist deshalb unzulässig. 


‘Zur Berechnung des Schallfeldes vorwärts bewegter Propeller geht man aus von GI. (24) 
für z schraubenlinig bewegte Quellen. Der Fortschrittsgrad der Quellen Damen wie oben 


= In = - ‚ die Absolutgeschwindigkeit der Anströomung V’=YV} + R?o?. Um mittels 


een mit dem Propellerelement mitbewegter Singularitäten die Randbedingungen am 


Propellerblatt befriedigen zu können, müssen zunächst in Gl. (24) zwei Dipolbildungen vorge- 
nommen werden. Zur Darstellung der Dickenwirkung des Profils benötigt man einen Dipol des 
Geschwindigkeitspotentials, dessen Achse stets in die Richtung y der schraubenförmigen Be- 
wegung zeigt.. Zur Darstellung der Anstellwirkung benötigt man andererseits das Feld eines 
Druckdipols, dessen Achse stets senkrecht zur momentanen Bewegungsrichtung steht. Die in 
Richtung der momentanen effektiven Anströmung zeigende Koordinate werde mit t,, die senk- 
recht dazu stehende mit n, bezeichnet. Man erhält aus Gl. (24) für den Dipol x- (Dipolachse 
in Bewegungsrichtung) bei Vernachlässigung quadratisch mit dem Abstand abnehmender Glieder 


4- ze dh, = (3022 2 4222) a, = (2Hesiny eur, - 


Om dt, 9001, ö öö R 
O,dtiz2o tg 9 ) . ( mzMr 
ee ee a A LER PR de 
2rr V 2 cos © Yig? O+ a VtgO +x? = 2 n YgO+r: (34). 
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n°c {0} 


-exp imzo (1 + 


Ebenso erhält man aus Gl. (24) für den Dipol y, (Dipolachse senkrecht zur Bewegungsrichtung) 
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NE RER ee 
DE m2 Mn exp |imzo(t en Vera N] 
Vg?O + 0+x #°c 7) 
Quadratisch kleine Glieder sind wieder vernachlässigt. 
20* 


Stellt wiederum y_ das Geschwindigkeitspotential und x, die Druckverteilung dar, so 
findet man zur Befriedigung der Randbedingungen am Propeller wie in (25) und 28) 
€ di,= DI V und Dam [4 1 5 ” * ” ” ” » ” ” [_! “ “ (36). N 
. . “ . . : * . t V ist 2 
Bei gleichförmig gerader Mitbewegung des Aufpunktes mit der Geschwindigkei z ist gemäß 
der een Eulerschen Gleichung der Zusammenhang zwischen Schalldruck und Ge- 
schwindigkeitspotential gegeben durch ERDE 


Del er 87. 


Wenn man mittels Gl. (37) in Gl. (34) zum Druck ‚übergeht und über den Propellerradius inte- 
griert, so erhält man analog zu der beim Propeller im Stand durchgeführten Rechnung im Fern- 
feld des dicken angestellten vorwärts bewegten Propellers für den Schalldruck 


M(z - ©)- Ms —n,)?+ x: +9). 
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Gl. (38) stellt mithin die umfassende Gleichung für das Schallfernfeld eines Propellers dar. Das 
Integral muß im allgemeinen Fall numerisch ausgewertet werden. 


d)Der angestellte Propeller im Standlauf (2. Näherung) 


Bei der Berechnung des Schallfeldes eines im Stand laufenden Propellers dürfen die mit 
der Strahlgeschwindigkeit nach hinten abschwimmenden Wirbel — wie oben schon gezeigt wurde 
— nicht mehr vernachlässigt werden. Man ersetzt deshalb in diesem Fall den Propeller durch 
einen mit der Strahlgeschwindigkeit V„ von vorn 'angeblasenen Propeller und berechnet das 
Schallfeld in einem relativ zum Propeller ruhenden Bezugssystem nach den Formeln (34), (35) 
und (38) für den vorwärts bewegten Propeller. Die Zusatzgeschwindigkeit V„ werde über den 
Strahlquerschnitt als konstant angenommen. Sie berechnet sich aus dem Propellerschub 75 
in hinreichender Näherung zu 
SER 


ae RE ER: (39). 


on Rz 
Die Geschwindigkeiten der gebundenen und freien Wirbel werden damit richtig wiedergegeben, 
iund das Schallfeld kann im Inneren und auf dem Rand des Schraubenstrahls in guter Näherung 
errechnet werden. Da die Luft im Außenraum ruht, müßte man mit den so auf dem Strahlrand 
errechneten Drucken als Randwerten mittels der Kirchhoffschen Formel oder der Singu- 
laritätenmethode das Druckfeld im ganzen Raum bestimmen. Dieser in der praktischen Durch- 
führung sehr komplizierte Weg soll hier nicht beschritten werden. Es wird vielmehr versuchs- 
weise das Feld des mit der Geschwindigkeit Y, angeströmten Propellers im ganzen Raum als 
Näherungslösung für den angestellten Propeller im Standlauf genommen. 

Man rechnet damit so, als ob auch außerhalb des Propellerstrahles übrall die Geschwindig- 
keit V, herrscht. Diese fehlerhafte Annahme führt bei der Phasenbestimmung zu entscheiden- 
den Fehlern. Bei der Berechnung des effektiven Schalldruckes dürfte die Auswirkung dieses 
Fehlers aber nicht so wesentlich sein, solange Y, < ce bleibt. Diese Bedingung ist im allgemeinen 
erfüllt. Das Schallfeld des angestellten Propellers im Standlauf wird damit ebenfalls durch 


; ; V, i 

Gl. (38) wiedergegeben, wobei M = == aus Gl. (39) bestimmt werden muß. Daß die so erhaltene 
Näherungslösung mit der Wirklichkeit recht gut übereinstimmt, zeigt der im folgenden durch- 
geführte Vergleich mit Experimenten. 


4. Vergleich zwischen Theorie und Experiment 


Zur Prüfung der Theorie liegen auswertbare Messungen an Flugzeugen nicht vor. Dagegen 
wurden von Ernsthausen[3] an Modellschrauben im Standlauf sehr sorgfältige experi- 
mentelle Messungen durchgeführt. Verwandt wurden für die Messungen zweiflügelige Modell- 


schrauben von 30. cm Durchmesser mit symmetrischen NACA- 


‘wurde für den unangestellten Propeller nach Gl. (26 


durchgeführt. Die Ergebnisse sind in Abb. 1 bis 3 dar- 
‚ gestellt. Da die Messungen in einer Entfernung von 
2 m durchgeführt wurden, mußte für den angestellten 
Propeller zur Berechnung des Schallfeldes das in (35) 


werden. 


Abb. 1 zeigt die für die Propellerebene gemessenen 
‚und berechneten effektiven Schalldrucke eines nicht 


Bild.1. Maximaler Schalldruck einer nicht- 
angestellten Schraube (Profil NACA 630012)’ 
in Abhängigkeit von der Umfangsgeschwin- 


digkeit 
Bild 2. Richtcharakteristik einer nichtange- =————ı=S1 
stellten Schraube im Standlauf (Profil NAOA — — —m=-2 | Meßwerte 
630012) für die ersten vier Hrmonischn re, — — — ——— m=3 ( Ernsthausen 
r (MR, = 0,8) en 
—m=1) oben: gerechnet EEE e 
REN Sa ggwten:Miwere —. 0 [--- m= (3) gerechnet 
Ernsthausen — m m = (4) 
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Bild 3. Schalldruckverteilung des Grundtones einer angestellten Schraube im Standlauf für verschiedene, Umfangs- 
geschwindigkeiten (Profil NACA 630006) 


gerechnet 


— nu u SOMESSEN 


angestellten Propellers für die ersten 4 Harmonischen in Abhängigkeit von der Umfangs- 
geschwindigkeit. Schon bei Machschen Zahlen von 0,6 liegen die Rechenwerte fast um 
den Faktor 2 zu hoch. Die Abweichungen dürften dadurch begründet sein, daß bei dem im 
Standlauf drehenden nicht angestellten Propeller die effektive Anströmgeschwindigkeit am Profil 
infolge Mitdrehens der Luft im Schraubenkreis sehr stark herabgesetzt wird. Die Größe der 
Quellkonstante ( müßte deshalb entsprechend verkleinert werden. Dagegen ist beim angestellten 


sehrau Profilen mit Dickenverhältnis 
Er ‚bis 18% und Anstellwinkeln bis 12°. Gemessen wurde die Schallfeldverteilung See 
_ ersten 8 Harmonischen bei Umfangsgeschwindigkeiten von 200 bis 400 m/sec (0,6 < Mps <1 Be 
Soweit die Messungen Er nsthausens in [3] veröffentlicht sind, wurden sie mit der 
- hier vorliegenden Theorie verglichen. Die Rechnung 


und für den angestellten Propeller nach Gl. (38) und (39) - € 


und (38) vernachlässigte Nahfeldglied mit berücksichtigt 
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Merbt u. Billing, Der Propeller als roti | | 
Propeller im Standlauf ein Mitdrehen der Luft in geringerem Maße zu erwarten. 


bewegten Propeller ist die Änderung der effektiven Anströmung im Rahmen unserer li 1ear 
sierten Theorie ganz zu vernachlässigen. = 


- Bi.‘ 
’ - ne 3 £ 
Abb. 2 zeigt die berechnete und die gemessene Richtcharakteristik der gleichen Modell- 


Schallfeldverteilung wird in den einzelnen Harmonischen richtig wiedergegeben. 


Abb. 3 zeigt die Richtcharakteristik für den Grundton eines angestellten Propellers im 
"Standlauf bei verschiedenen Umlaufgeschwindigkeiten. Auch hier liegen wegen des Mitdrehens { 
- der Luft im Schraubenkreis die Rechenwerte noch etwas über den Meßwerten, doch ist die Über- 


einstimmung wesentlich besser als im Fall des unangestellten Propellers. Die räumliche Schall- 


feldverteilung wird für Mr, = 0,6 und Mr, = 0,7 gut wiedergegeben, dagegen treten bei 


Mr, = 0,9 Abweichungen auf. Diese starken Abweichungen beihohen Mach schen Zahlen er- 


klären sich daraus, daß zur Druckbestimmung am Propellerprofil o„-Werte benutzt werden 
mußten, die bei kleinen M a c h- Zahlen gemessen waren. Durch Berücksichtigung derPrandtl- 


\ + C 
Kontraktion wurde daraus ein für die Rechnung gültiges ag = . ya gewonnen, was nach 
j V EI 
c®? 


den bisherigen Erfahrungen im Bereich hoher Unterschallgeschwindigkeiten nur bis etwa 


Mr, = 0,75 zulässig ist. Darüber hinaus erhält man (vgl. z.B. Helmbold [10]) zu große 


- Ca-Werte. 


Hinsichtlich der Obertonverteilung ist für den angestellten Propeller ein Vergleich mit 
Experimenten leider nicht möglich, da Meßwerte nur für Mr; = 0,9 vorliegen. 

Für den Fall des vorwärts bewegten Propellers können die vom Propeller induzierten Zu» 
satzgeschwindigkeiten (radiale und axiale Zusatzgeschwindigkeit der Luft im Propellerkreis) 
vernachlässigt werden. Die Ergebnisse der Theorie (nach Gl. (38)) dürften deshalb mit der 
Praxis wesentlich besser als im Standlauf übereinstimmen. Als Beispiel wurde das Schallfeld 
eines Flugzeugpropellers natürlicher Größe berechnet 1), Die Richtcharakteristik wurde getrennt 
nach Dickenanteil (Abb. 4) und Schubanteil (Abb. 5) aufgezeichnet. Abb. 6 gibt die durch Super- 
position gewonnene gesamte Richtcharakteristik der angestellten vorwärts bewegten Schraube. 
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Bild 4. Schalldruckvortei- 


‚I 
de TR 


schraube für die ersten 4 Harmonischen für eine hohe Umfangsgeschwindigkeit. Die räumliche 


Ei 


Bild 5. Schalldruckverteilung 


lung des Grundtones einer des Grundtoneseiner vorwärts- 


vorwärtsbewegten Normal- 


bewegten Normalluftschraube 
en. Dicken- — Schubantell 
8) Schraube (% = 0,8; 2}. Eolranbe + a es 
MR = 0,8) Rs "rise 
b) Schraube II (A = 0,5; b) DanreRDe Ne: 2 
Mrs = 0,8) Rs ’ 


Bild 6. Schalldruckverteilung des 
Grundtones einer vorwärtsbeweg- 
ten Normalluftschraube — Ge- 
samtschalldruck 
a) Schraube I (4 = 0,2; 
Mrs = 0,8) 
b) Schraube II (A = 0,5; 


Mr; = 08) 


Entsprechend dem Überwie 


gen der Anstellwirkung ergibt sich eine Verstärkung des Schall- 
druckes n 


ach rückwärts. Das der Dickenwirkung des Propellers entsprechende sekundäre Maxi- 
mum in Flngrichtung (bei © = 30°) tritt hier nur sehr schwach auf, da für die Normalluft- 
ı) Entwurf einer der Bestschraube angenäherten Normalluftschraub 


i e nach Helmbold ([10], C, 2) für 
4 == 0,2 und A == 0,5. Mach-Zahl der Propellerspitzen Mr s = 08, Leistung 1000 PS. 
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schraube ein dünnes Profil gewählt wurde. Dickere Profile ergeben — unter sonst gleichen Be- 
triebsbedingungen — eine Vergrößerung der Werte in Abb. 4 und entsprechend eine Verschärfung 
des Nebenmaximums in Abb. 6. ne 

Die Maximalwerte des Schalldrucks entsprechen mit pr = 1400 ub : m Schalldrucken, wie 
sie im Flugbetrieb bei Propellern mit hohen Umfangsgeschwindigkeiten beobachtet werden. Für 
kleine Umfangsgeschwindigkeiten (etwa M. Rg = 0,6) ergibt sich bei gleicher Motorenleistung 
ein wesentlich niedrigerer Maximalwert von 900 ud. m. | 

Allgemein läßt sich aus der Formel (38) noch folgendes ableiten: ' 

Propeller ohne Anstellung ergeben zweiblättrige Richtcharakteristiken mit entsprechend 
dem Fortschrittsgrad in Flugrichtung vorgeneigtem Maximum (Abb. 4). Dagegen gibt der 
zweite Term der Gl. (38) vierblättrige Charakteristiken, das Maximum des abgestrahlten Schalls 
liegt hinter die Schraubenebene verschoben (Abb. 5). Beide Anteile sind etwa von gleicher 
. Größenordnung. Mit steigender Blatt- bzw. Obertonzahl nimmt der Einfluß der Anstellung im 
Vergleich zur Dickenwirkung ab. 
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Zur Theorie der tragenden Linie mit periodischer Zirkulation 
Von Herbert Bilharz in Freiburg i. Br. 

. Es wird die zweite Grundaufgabe der Prandilschen Theorie der tragenden Linie für eine vorgegebene 
periodische Anstellwinkelverteilung gelöst und die erhaltenen Ergebnisse an zwei Beispielen veranschaulicht. 


The second fundamental problem of ihe Pranditl theory of the. trailing line is solved for a given 
periodical distribution of the angel of incidence, and the resulis obtained are illustrated by two examples. 


Le .second problöme fondamental de la ih£orie de la ligne portante selon Prandtl est rösolu pour 
une distribution periodique donnde de l’angle d’incidence. Les resultats trouvds sont illustres par deux 
examples. 

Pemaerca BTOpaA OCHOBHaA 3anaya Teopun INIOCKOrO Kpuna Ilpanırıa ma sanaunoro 
mepuonuyeckoro pacnpeneneunun yrıoB atakı. llonyuenmse pesyıbTaTH WINOCTPHPyITcH 
Ha AByX IIpuMepax. 

1. Einleitung 

Als L. Prandtl!) vor 31 Jahren in seiner ersten Mitteilung zur Tragflügeltheorie die 
tragende Linie mit periodischer Zirkulation als einfachstes, seiner Methode zugängliches Bei- 
spiel behandelte, schien jenes Problem nur theoretisches Interesse zu besitzen. Inzwischen hal 
es längst in mehrfacher Hinsicht praktische Bedeutung erlangt, worüber u. a. eine Patentschrift ?) 
aus dem Jahre 1940 unterrichtet, die eine periodische Verformung der Hinterkante von um- 
strömten Körpern zum Gegenstand hat, um größere Wirbelballen in einzelne kleinere Wirbel- 
zöpfe aufzulösen. Es schien daher angebracht, jenes Beispiel Prandtls, zumal es a.a. 0. 
nur kurze Erwähnung findet, wieder aufzugreifen und in seinen Einzelheiten weiter zu verfolgen, 
insbesondere den induzierten Widerstand und das Abwindfeld zu bestimmen. Dies ist das Ziel 
der vorliegenden Mitteilung. Zwei Beispiele werden zur Belebung des dargebotenen Formel- 
apparates explizite durchgerechnet. 


2. Bestimmung der Zirkulation 
Im unbegrenzten Raume sei bezüglich eines kartesischen (z, y, 2)-Koordinatensystems eine 
stationäre, zur (x, 2)-Ebene parallele, reibungsfreie und inkompressible Strömung gegeben. Der 
ı)L.Prandtl, Tragflügeltheorie, 1. Mitteilung. Nachr. Ges. d. Wiss., Göttingen, Math.-phys. Kl. 


1918; S. 451—477; insbes. S. 471—473. 
2) Patentschrift Nr. 695428, Kl. 74d, Gr. 6,, (Berlin 1940). 
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Geschwindigkeitsvektor habe im Unendlichen den Betrag V; die Dichte des Mediums sei e = 
const >0. Die y-Achse trage eine kontinuierliche Wirbelbelegung y (y), die eine. vorgegebene 
örtliche geometrische Anstellwinkelverteilung & (y) hervorbringe, welche mod 22 periodisch 
(0 <1 < ©) und durch die konvergente Fouriersche Reihe 


&(y) -3+ (mon + busin.n 2) Re: 5 2 AR EEE 


dargestellt sei. En } ! 
Unter den, der Prandtlschen linearisierten Theorie zugrunde liegenden Voraus- 
setzungen ®) folgt die zu bestimmende Zirkulation y(y) aus der Integrodifferentialgleichung - 


1a re, (22 
Ay) gzy Fre ga nV PER TE ER e ( ) 


Eu 


mit dem Cauchyschen Hauptwert an der Stelle y = y i j “ 
Zur eindeutigen Auflösung dieser Gleichung nach y(y) wählen wir den Reihenansatz: 


4, x nY ey 
y(y) -av|5 +2 (Ancos NT + B„sin „“r REN (2.3).- 
Mit der Substitution n= y— y’ und Beachtung von 
elanan, PR Man=0... ren: (2.4) 


für u >0, folgt nach einiger Zwischenrechnung ‘) 


1 ddr dy _MVS ( ny any : 
TE TE A„cosn 7 + B„sinn RENTE (2.5) 


und nach (2.1) und (2.2) für die Koeffizienten von y in (2.3) 


b 
VATER DER m Bn= —i— WER eo 
N) 4 


3. Auftrieb und Widerstand 


Mit einer fest gewählten Längeneinheit als Bezugstiefe hat der Auftrieb (z-Komponente 
der resultierenden Kraft) über eine volle Periode den Betrag 


[i 


Aey= 2 P:21=0 | yiy)ay re. (3.1); 
-ı 
unabhängig von der Periodenlänge ergibt sich hieraus ein dem mittleren geometrischen Anstell- 
winkel a,/2 proportionaler Auftriebsbeiwert 


Der Betrag des induzierten Widerstandes (x-Komponente der resultierenden Kraft) über 
ein volles Periodenintervall ist 


2 
wen=ct ve 21= 0 |yimoray er REN ER (3.3), 


1 
worin die Zirkulation y aus (1.3) und der Betrag 


®) L. Prandtl, loc. cit. 8. 466-467. 

*) Hierbei vorzunehmende Vertauschung in der Reihenfolge von Integration und Summation ist nach 
einem Satze von Ü. Arzelä gestattet, falls wir voraussetzen, daß die Folge der Partialsummen des Inte- 
granden, die gegen eine intergrierbare Funktion konvergiere, selbst integrierbar und gleichmäßig beschränkt 
sei. Vgl. etwa F. Hausdorff, Beweis eines Satzes von Arzelä, Math. Ztschr. 26 (1926), 8. 135—137. 


S Zuhilfenahme von (2.5) een an Y legen der Orthogonalitätscigenschaften 4 der -trigono- R. a 
nn ae Euoksionen erhalten wir das. ne ee | ! A A 


\ w N fi , [= Se ee Ce TE Eon 


. N einen von ne Periode RR Widrtandbeinert. Da Z, (a ar Z N) konvergiert, ist 
ie die le Reihe in 8. 5) a Bacon: konyergeut., ; : ER 


E } 4. Abwindfeld® en PR ed 
Rear Nach dem Geketz) von Biot und Sa vart ist allgemein der Betrag. w "des von einer E 
ER . Wirbelinie, s im m Aufpünkte iz ee een w en Bes me. 


. gegeben, worin 7 (s) die Zirkulation, ds das Tinlenelemeht: von sr den: Abstand ee ne A 
B.... _ vonds undg(s) den Winkel zwischen der Richtung des Wirbelelementes und seiner Verbindungs- Ko 
geraden mit P bezeichnet. Der Vektor tv steht normal zur Ebene ‚durch 7: und ur seine Orien- R 
F tierung um s stimmt mit derjenigen. von y(s) überein. IR 
4 Um hiernach das Abwindfeld zu berechnen, Bee wir zunsehet, mit er Bestimmung > 
. des Abwindes %y, der von der y-Achse (gebundener Wirbel) in einem Aufpunkte P(2,9;0) mit 
#: . &.>0 induziert wird. Hierbei kann die y-Koordinate des Aufpunktes Ra der Periodizität 
auf das Grundintervall —I <% <sI! beschränkt werden. 2 
Nach Einführung der Polarkoordinaten 


r=eyRr+y— N >0, y=arcigYt je v<d).. | (2) n 
folgt aus (4:1) mit beschränkt EN r+(y) aus (2. 3) und AU Integration 
u, (4,9) = ae 
j Erd RR I RS RE ER 
e Zar (& |eos In? (ätep+ 3] cospdy-+Bn j sin In# (Ztgy+ D) a) (4.3): 
$ n=1 . \ ir ; 
Re b3 EEE Er | | 
Br .. und daraus mit der weiteren Variablentransformation VER De: 
2 Sint=tgy -osıza; —5 sy <=) 3 (4.4) > 2 
schließlich | } 2 
uy (8,4) = 7; +2 (An cosn"F 4 Basinn® “ —dr Wo 
% 4o 2 . N &ofr. 


Da die modifizierte Besselsche Funktion zweiter Art und p-ter Ordnung, Ken 
welche mit Hilfe der modifizierten Funktion erster Art: 


4 wi I: So 1 p-+23 
ER m FA {RE — et 3 Be an ee (4.6) 


s=0 
für nicht negatives ganzzahliges p durch 


gr vi 
ri -5Q), 


A ee a (4.7) 


2c0sprn 
‚definiert ist, der Bassetschen Integraldarstellung °) 


) 
K ln (e) cos (£ Sin Tr) 
NS Cop? r 
0 
5) G.N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge 1922), S. 172; Funk- 
tionswerte daselbst 8. 698—713 und S. 737. 
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mit der Eulerschen Gammafunktion 


r(p+ Bin (>0), r(&) =Ya 


o 
TE 


I 
win_ 1A, Sr ("*#)] =y = . . (4.10). 
ar Din , K,ln ) Ay ecosn ) + B, sin n 1 (4.10) 
Für 20 ist 0O<CK,(£) <1, denn £K,(£) ist monoton fallend und 
Kl _ 
nn —t 
3E e 
bei festem % ist daher w, aus (4.10) von der Ordnung &"! und kann deshalb für das Folgende 
unberücksichtigt bleiben. a Y 
Die induzierte Geschwindigkeit v; (#, 5) der freien Wirbel bestimmt sich ebenfalls aus (4.1), 


indem wir zunächst parallel zur x-Achse, sodann über die y-Richtung integrieren und an der 
Geraden y= % den Cauchyschen Hauptwert nehmen: 


genügt, erhalten wir für p=1 und {=n—- aus (4.5) und (4.8) 


limöK,($)=1, lim 
t>1 (>> 


IN 1.LE) z dy 
v(% BT —I(1 nt ee I a BR I TE Te a (4.12). 
er rer) | 
Der erste Summand dieses Integrals ist aus (3.4) bekannt und hat nach (2.5) den Wert 
"V< rn =) 
= — > — Er A nA .15); 
v 41 Di n(Ancosn ; + B„sinn ; (4.13) 


es ist danach hinreichend, das Integral 


ER an ; dy 
v(%; — Mar eugi £ nr en a (4.14) 
NT ER F Yard 
auszuwerten, welches für &— » ebenfalls gegen v konvergiert. 
Mit (4.2) geht (4.14) über in 


© 3 o 
vi(z, 9-7 n(Ancosn + Basinn") | sin [Per] . . (4.15). 
Wal 
-3 


Wird zur Abkürzung für y 20 


5 dy 
Hy= | sin (xIgy) Ar 
-3 
a so folgt, da die zu (4.8) führenden Überlegungen auch im Limes p = 0 gelten, daß das 
ntegral 


Fi 
IdH 1 d 


2dr a (xtgy) cosy — j cos(F Sir DET (4.17) 
-# Ö 


existiert und gleich der modifizierten Besselschen Funklion K,(x) von zweiter Art und 
nullter Ordnung ist; wegen H (0) = 0 ist daher 


und es gilt nach Heaviside®) 


lim Hg=2 | Kddi=n. . (4.19). 
z> > Ö 


%)G. N. Watson, loc. cit. 8.388. 


Bild 1 Sam den Verlauf von u @ für En 
a EAU Dnu( a 
-—— — — en Wil 


‘und (4.20) den Abwind der freien Wirbel in Punkten der (z, y)- -Ebene: 
RE < u. .ny E 
y-Y v+%)= 212 n (& + un ”)) (Ancosn "Bari 7) .. (4.21). 


.schwungenen bzw. gekröpften) Profilstäbe 


Aus (4.15) und (4. 16), falls darin. en r 
gesetzt wird, erhalten wir sodann rc 


(420) 


[Area N +B, it 


und schließlich durch Zusammenfassen von (4.13) rk 


5. Beispiele 
- Als Beispiele für ein mod.2] periodisches &(y) wählen wir . 
a) eine reine Kosinusverteilung um den Mittelwert #0: 


a(y= tote (0, ON er . (5.a), 
b) eine gerade, EDIeLEBNEIRE, konstante Funktion: 
03 0o<|yl<y2, 
Cg+ 03 5 Sr AN Mi 
a7 für Sl N Ee 


wobei ohne Beschränkung der Allgemein- 
heit die Konstante c, größer als die Kon- 
stante c, angenommen werden kann. 
Diese Verteilungen geben die einfach- 
sten Modellvorstellungen für die in Bild.2 
dargestellten periodisch verformten (ge- 


von unendlicher Länge und konstanter’Tiefe. 
Die Fourier-Konstanten der 

Reihe (2.1) sind 

al id a 

net lürie ned, Isa 

beziehungsweise 


= _&a+% Pr _ 28 — 6) 
2 A a PATE 
ER —2(0,— 65) (5.1b). 
EEE EETT ER 
aaın41) = dn41ı = 0 für n=(, 1, ... Bild 2 
Durch die Gleichungen 
‘=6y Er (3.28) 
beziehungsweise durch 
EN i 
a) o= a (9.2b) 


sei ein mittlerer geometrischer Anstellwinkel & und ein Differenzwinkel 6, der das Maximum 
von |& — «| angibt, erklärt. Mit diesen Bezeichnungen und den aus (2.6) errechneten Koeffi- 
zienten lauten dann nach (2.3) die Se 


TyY 
vy(y)=nVla-+ "cos — (5.3a) 
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und 


cos an +1) 3 cos (4n +3)? 
(Y=rV SEN Ä i RE 
; ) \an+ 1 (1444177) an+ 31 +4n+3%) 


Für die Belegung (5a) ist die Zirkulation (5.3) längs der y-Koordinate ebenfalls eine um 
ihren Mittelwert auherde Kosinuslinie, deren Amplitude mit } monoton wächst. 


Im Falle (5b) ist mit «= ö, also og 0, die Zirkulation (5.3b) im Grundintervall 


(5:3b). 


—1<y <lfür verschiedene Werte von lin Bild 3 aufgetragen. Die Amplitude nimmt 'eben-- 


falls mit ! monoton zu; die Maximalwerte beider Verteilungen stimmen näherungsweise überein 
und in beiden Fällen findet ein Ausgleich der Strömung nach den Periodenenden hin statt; mit 
abnehmender Periodenlänge strebt der y proportionale effektive Anstellwinkel gegen den Mittel- 
wert a. 

Der Auftriebsbeiwert ist nach (3.2) 


Dur in — ER a: (5.58) 
( +2) 
und 
= air — a (5.5 b) 
= >n-+ (1 +(2n+ y%) 


also unabhängig von & und quadratisch in 6. 


Bild 3 


Bild 4 
Die Abhängigkeit von list für (5.5a) leicht festzustellen, denn diese Funktion hat an der 
2 
Stelle = T ein endliches, positives Maximum: max c„= eZ und strebt, wie es sein muß, 


nach Null bei den Grenzübergängen Z>0 und Z>o. Ein ähnliches Verhalten zeigt die un- 
endliche Reihe in (5.5b). Mit Hilfe der Eulerschen Summenformel”) erhält man bis auf 


einen von l unabhängigen Fehler, der dem absoluten Betrage nach kleiner als 0,0017 gemacht 
werden kann, 


oo 


el 10T RE 
on li+entnZ) 145) s(1+37) (1437) 


1 41 1 1 24m? 
77 ar \ ar 


% AR Vgl. K. Kno’pp, Theorie und Anwendung der ‘unendlichen Reihen, 3. Aufl. (Berlin 1931), S. 540 
is 3 


(5.6. 


{ 


(ya Le he ers 


4s wel 1,b) a als NE von U graphisch darges Ute: 
B en konstanten geometrischen Anstellwinkel Ms) b) ist w ungefähr doppelt 


wie im Falle & a); das a liegt näherungsweise bei i=3> 7 . (Im Falle (5b) konnte 


durch Experimente, bei denen I auf das Intervall 0,1 <I < 1, 6 beschränkt lol das Anwachsen | 


des induzierten Widerstandes qualitativ festgestellt werden. )v 


Aus (4. 10) folgen die von den gebtindenen. Wirbeln induzierten Geschwindigkeiten. 


RE Ban 
RD 2 2 29 
| Ita) ER en 


(5. a 


RR 


Hieraus ergeben sich für = 0, a= öund Berkchene Werte von = -i= - const. >0 Kurven- 
scharen über &, die monoton fallend zwischen —- % und ! - verlaufen; der asymptotische Wert 
33 ist bei den Abszissen. &=2], also Rrans eine Periodenlänge von der yrAchse entfernt, 
nahezu erreicht. RE 
Mit Hilfe von (4.21) erhalten wir schließ- 
lich die von den freien Wirbeln en 


Ge schwindigkeiten 


Y; IT £ 
7-87, 5 = (5.82) 
Ber an 

und 
3-22. E Tun 

| Han), 

 a+Hlen EZ n 
| an. (4n +) 
1+(4n +37; 
(5.8b) 


Für 7=0Osind in Bild 5 die Funktionen v,/V aus (5.8a) über und dem Parameter 2 wieder- 
gegeben. Der Abwind ist von & unabhängig, proportional mit öund umgekehrt proportional zu |. 
Wie die Abbildung erkennen läßt, wird der im Unendlichen herrschende Strömungszustand 
bei kurzen Periodenlängen rascher angestrebt als bei langen und ist etwa für =], also un- 
gefähr eine halbe Periodenlänge von der y-Achse entfernt, nahezu erreicht. Das gleiche gilt 


für (5.8b). 
Eingegangen: 24. März 1949. 


(tet) eos(An+ 1’ a Ze 


Prof, Albert Rohrberg (Abteilungsleiter a. d Tue: 
nieurschule Gaußin Berlin), Graphische Funk- 
tionstafeln, enthaltend die Mantissen der de- 
kadischen Logarithmen und die goniometrischen Funk- 
tionen von Minute zu Minute auf vier Stellen. 30 8. 


"Berlin 1949. Fachverlag Schiele und Schön. Preis 


brosch, 4,50 DM. 


Die graphischen Tafeln enthalten Doppelskalen, in 


deren eıner die Nummeri bzw. Winkel, in deren anderer 
nis i nten ‘der Winkel. 


die Logarithmen, Sinus bzw. Tange. 2 
aufgezeichnet sind. Dabei liegt immer einem gleich- 
mäßigen Skalenteil ein leichmäßiger gegenüber, 
aber keine der Skalen ist ihrer ganzen Länge nach 
gleichmäßig; sondern, je nachdem sich die günstigste 
Ablesung ergibt, teils gleichmäßig, teils ungleichmäßig. 
Die Skalen sind natürlich gebrochen, z. B. die für die 
Logarithmen achtzigmal. Die et er dieser 
einundachtzig Teilskalen erfordert acht Seiten. Die 
Maßeinheiten sind so gewählt, daß man ohne Inter- 
polation vier Stellen, mit Interpolation fast überall 
5 Stellen ablesen kann. Aber der Zweck dieser Tafeln 
soll es gerade sein, jede Interpolation zu ersparen. Da- 
für wird der große Vorteil gewöhnlicher vierstelliger 
Tafeln, daß es möglich ist, sämtliche Logarithmen 
auf zwei benachbarteri Seiten unterzubringen, so daß 
bei ihrem Gebrauch ein Umblättern unnötig ist, auf- 
gegeben. Dem gegenüber scheint mir die Arbeit, die 
das — vor allem bei Benutzung von Proportionaltäfel- 
chen — leicht im Kopf auszuführende Interpolieren 
macht, ein geringer Nachteil, zumal auch. das Auf- 
suchen und Ablesen einer Zahl in einer Skala, selbst 
wenn sie so gut und übersichtlich ist, wie in den vor- 
liegenden Tafeln, umständlicher ist, als das einer ge- 
druckten Zahl in einer Tabelle. 


Dresden. 


Dr. W. Volkmann (Oberstudienrat a.D.), Ele- 
mente physikalischer Experimentier- 
kunst. 2. Auflage. VIII + 173 S.m. 141 Abb. Bonn 
1949, Ferd. Dümmlers Verlag. Unveränderter Neu- 
druck der Auflage von 1939. Preis brosch. 5,80 DM. 

Das vorliegende Buch beschreibt einfache Versuche, 
die im Physikunterricht der Schule ausgeführt werden 
können und gibt dabei ausführliche allgemeine Hin- 
weise für eine zweckmäßige Vorführtechnik. Es wird 
allen denen eine sehr erwünschte Hilfe bedeuten, die 
sich die Grundelemente physikalischer Unterrichts- 
versuche aneignen wollen oder die eine Unterrichts- 
sammlung neu aufbauen müssen, 


Dresden. 


Willers. 


Recoknagel, 


Dr.-Ing, habil. K. Hirschfeld (o. Prof. a. d. Techn. 
Hochschule Aachen, Die Temperaturver- 
teilung im Beton. 154 S. mit 173 Abb, und 
15 Zahlentafeln. Berlin, Göttingen, Heidelberg 1948. 
Springer-Verlag. Preis brosch. 36,00 DM. 

Die Frage nach der Temperaturverteilung in Beton- 
körpern wird erstmalig in einer physikalisch und mathe- 
matisch einwandfreien Form angeschnitten und für be- 
stimmte Fälle gelöst, Ausgang der Untersuchung ist 
die allgemeine MRraaleihungschelchung, die unter Ein- 


führung der Temperaturleitzahl « = n und des La- 
W 


6 

angeschrieben wird, Durch Beschränkung auf Han 
parallele Betonplatten und Betonzylinder von unend- 
licher Ausdehnung und durch Idealisierung der Pro- 
blemstellung wird das dreidimensionale Problem in hin- 
reichender Annäherung auf ein eindimensionales zu- 
rückgeführt und somit überhaupt eine numerische 
Lösung ermöglicht. Da für die praktische Anwendung 
der Rechenaufwand immer noch sehr umfangreich ist, 


placeschen Operators in der Form = =a49+ 


Rücksicht auf zahlenmäßige Berechnungen 


‚on F EEE Is Par er i 
k : 'Buchbesprechun gen fr hi j 3 er nr Ele 


BUCHBESPRECHUNGEN 


hat der Verfasser die Ergebnisse in zahlreichen Schau- 
bildern 


Aufgaben Bee wird. i 
Dresden, 


Dr. W Gröbner und Dr. Nikolaus Hofreiter, 
Integraltafel, Teil I: Unbestimmte 
Integrale. VIII+ 166 8. Wien und FU 
1949. Springer-Verlag. Preis brosch. 18,00 DM. 


Henn, 


Die vor ; Integraltafel stellt eine Neubear- 
beitung an krwei 


terung des Notdruckes dar, den die 


Verfasser im Jahre 1944 herausgebracht haben. Es ist 


anzuerkennen, daß sie sich der großen Mühe dieser 


Neubearbeitung unterzogen 

der mathematischen Literatur Deutschlands ausfüllt; 
denn seit 1900 sind meines Wissens in deutscher 8 

bisher nur zwei Integraltafeln erschienen, zwar 
1906 die ‚‚Tafeln unbestimmter “von G. Petit 
Bois und 1935 die ‚‚Integralformeln für Ingenieure und 
Studierende‘‘ von Naske, die ja kaum größeren An- 


sprüchen n. Die vorliegende Integraltafel glie- 
pad sich Ba Abschnitte: I. Rationale Integranden; 
IL braische irrationale In n; Tran 
szendente Inte . Die Schreibweise der Integral- 
formeln wurde quent der funktionentheoretischen 
Auffassung paßt, doch haben die Verfasser darüber 
hinaus auch den Forderungen der Praxis Rechnung 
getragen, indem sie die Formeln, soweit es erwünscht 
schien, außer in der funktionentheoretischen Schreib- 
weise auch in einer für die numerischen Auswertung 
im Reellen geeigneteren Form angegeben haben. Zur 
Transformation der elliptischen Integrale auf die Le- 
gendresche Normalform benutzen die Nee mit 

e qua- 
dratischen Transformationen. phisch stellt die 
Integraltafel von Gröbner und Hofreiter ein Kuriosum 
und ein Meisterwerk dar: vom Vorwort und Inhalts- 
verzeichnis abgesehen ist: nichts in Lettern gesetzt, 
sondern alles handgeschrieben und durch Steindruck 
vervielfältigt. Der r wird an der im einzelnen wie 
gen: sauberen Ausführung seine Freude 

aben, | 


Freiberg Sa. 


Grüß. 


Hans von Mangoldt — Konrad Knopp, Einfüh- 
rungin die Höhere Mathematik. Bd.I 
685 S. mit 112 Abb., Bd. II 634 S, mit 108 Abb., Bd. 
III 618 S. mit 103 Abb. Stuttgart 1948. Verlag S. 
Hirzel. Preis pro Pand 24,80 DM. 

Unter den so mannigfaltigen, mehr oder weniger 
strengen, leichter oder schwerer zu verstehenden, mehr 
den Anwendungen oder mehr der Theorie zugewandten 
Lehrbüchern der Höheren Mathematik werden die drei 
Bände von. v. Mangoldt-Knopp immer ihren bevor- 
zugten Platz behalten; sie ve n ihn der glück- 
lichen Mischung von Strenge und Anschaulichkeit, von 
abstrakter Theorie und lebendiger Anwendung. Von 
der siebenten Auflage an sah der Verfasser keine Mög- 
lichkeit mehr, sein Werk zu verbessern, es hat also 
offensichtlich seinen asymptotischen Idealzustand 
schon soweit erreicht, daß höchstens noch etwaige 
Druckfehler auszumerzen bleiben. So ist auch die 
neunte Auflage ein fast unveränderter Abdruck der 
achten. 

Freiberg. Grüß, 

Dr. phil, habil. Bruno von- Freitag gen. Loringboft, 
Gedanken zur Philosophie der Ma- 
thematik. 53 S., Meisenheim am Glan 1948. West. 
kulturverlag Anton Hain. Preis kart. 3,20 DM. 

"Das lesenswerte und gedankenreiche Büchlein hier zu 
erwähnen — eine Besprechung verbietet sich aus räum- 
lichen Gründen —, erscheint vor allem deshalb gerecht- 


eh. a Zn 


NAT NER 


eo, MW a rn u / 500 er 
ma WILL, 


 Bd.20Nr.100%t.1940 


U 


Dr.-Ing. W. Meyer zur Capellen, 


Z,angew Math.Mech, 


. [RZ ), 
v „ ” 


k 


fertigt, weil'darin das Problem der An wendbar- 


‚keit der Mathematik auf die Wirk- 
lichkeit am ausführlichsten und in einer recht: les- 


baren Form behandelt wird. Die philosophische Grund- 
lage dazu bringt der erste Absohnitt, der die Seinsweise 


. der mathematischen Gegenstände erörtert. An diesem 
‘ bleibt freilich zu wünschen, daß der leicht irreführende 


Begriff des Fiktiven eindeutig festgelegt wird. Aus dem 
letzten die „menschliche Seite der Mathematik‘ be- 


handelnden Abschnitt sind vor allem die anregenden 
. Betrachtungen zum Methodenproblem hervorzuheben, 


Dresden. M. Draeger. 

Mathemati- 
sche Instrumente (Mathematik und ihre An- 
wendungen in Physik und Technik. Herausgegeben von 


Prof. Dr. E. Kamke und Prof. Dr. A. Kratzer, Bd.1). 


Dritte, ergänzte Aufl. X + 339 S. mit 250 Abb. Leip- 
zig 1949. Akademische Verlagsgesellschaft, Geest u. 
Portig Kom.-Ges. Preis geb. 23,00 DM. 


-Von der vorliegenden dritten Auflage des Buches, 


dessen erste in dieser Zeitschrift Bd. 21 (1941) 8.63 aus- 


‘ führlich besprochen wurde, waren zu Beginn des Jahres 


1945 sämtliche Bogen ausgedruckt, doch konnte das 
Werk in jener Zeit nicht erscheinen. Der Inhalt war 
auf den damaligen Stand gebracht; unter anderem wa- 


‚ren eine Reihe von Abbildungen, durch deren Güte sich 


das Buch von Anfang,an auszeichnete, durch andere, 
die neuere Konstruktionen wiedergaben, ersetzt. Nach 
1945 wurden dann zahlreiche neue, während des Krie- 
ges meist im Ausland konstruierte Apparate bekannt, 
vor allem die in den USA gebauten Rechenautomaten 
und die Neukonstruktion der Integriermaschine von 


Bush, auch konnte erst jetzt über die Entwicklung 


in Deutschland während des Krieges berichtet werden. 
Da der Verf. auf all dieses eingehen wollte, brachte er es 
in einem Nachtrag, in dem außerdem einige neuere har- 
monische Analysatoren und Synthesatoren behandelt 
werden. Auch das’ Schrifttumsverzeichnis, das jetzt 
346 Nummern umfaßt, wurde durch einen Nachtrag er- 
gänzt. Vielleicht sind in diesem Anhang die Rechen- 
‚automaten, bei denen nach dem Kriege eine äußerst 


” Eingegangene Bücher 


USA, sondern auch in England, Frankreich und der 
USSR, dabei etwas zu kurz gekommen. Das läßt sich 
aber bei einer Neuauflage, bei der auch dieser Nachtrag 


eingearbeitet wird, leicht nachholen. 
Jas während des Krieges ins 


Nachschlagewerk sein. Möge es 
mansichin Deutschland, das bis 1945 für die Konstruk- 
tion vieler Gruppen mathematischer Instrumente füh- 
rend war, wieder lebhafter als das in den letzten Jahren 


‚möglich war, an der Entwicklung auf diesen Gebieten 


beteiligt. i 
Dresden. Willers. 


Dr. WalterLietzmann (Prof. a.d. Univ. Göttingen), 
Elementare Kegelschnittlehre, eine 


‚Einführung in die Methoden der Geometrie. 171 S. mit 


104 Abb. Bonn 1949. Ferd. Dümmlers Verlag. Preis 
brosch. 11,80 DM. 


Ein inhaltsreiches Buch, in dem die Eigenschaften der 
Kegelschnitte aus den verschiedenen möglichen Defi- 
nitionen so abgeleitet werden, wie es jeweils am zweck- 
mäßigsten undeinfachsten ist. Die einprägsame Gesamt- 
übersicht und Kai rin einge wird dadurch et- 
was erschwert. Da den einzelnen Abschnitten Übungs- 
aufgaben beigegeben sind, dürfte die Herausstellung 
dieses Vergleichs eine der wesentlichen en sein, 
die der Leser nach der Durcharbeitung des Buches als 
einen Querschnitt unter anderem Gesichtswinkel zu 
lösen hätte. TR ! 

In den einzelnen Kapiteln werden die planimetrische, 
stereometrische, analytische, affine, perspektive, pro- 
jektive und gruppentheoretische Definition besprochen; 


der Äquivalenznachweis dieser Definitionen wird er- - 


bracht. — Wenn als Untertitel „Eine Einführung in die 
Methoden der Geometrie‘‘ genannt wird, so wäre nach 
dieser Richtung hin eine etwas Sag er ES pielle 
Darstellung— widsie am Schluß des ersten Äbschnittes 
begonnen wird — wünschenswert gewesen. 


Dresden. M. Draeger. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei.der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Dr. Walter -Lietzmann (Prof. a. d. Univ. Göbtingen). 
Das Wesen der Mathematik. (Sammlung 
„Die Wissenschaft‘‘, Band 102; Herausgeber Prof. Dr. 


Wilhelm Westphal.) VII/168. 8. mit 41 Abb. Braun- 
schweig 1949. Verlag Friedr. Vieweg & Sohn. Preis 
karton. 8— DM. 


Druckfehlerberichtigung 


zu L. Collatz, Differenzenverfahren zur 
numerischen Integration... Z. angew. 
Math. Mech. 29 (1949), S. 199-209. 


Im ganzen Artikel lege man ‚‚Differenzenschemaver- 
fahren‘‘ statt ‚‚„Differenzenverfahren‘‘. 


8.199, 2. Absatz, Zeile 7: „Ersetzens‘‘ statt „Ein- 
setzeng‘', 

Gleichung (7): „FD (&“ statt „OrV(H“. 

Gleichung (11): Rechts füge man hinzu „+ R,'“. 

Gleichung (13): Das Integral über u hat 0 (nicht «) 
als untere Grenze. 

Gleichung (14): In der Summe deroberen Zeile „‚f(2,)‘“ 
‚statt ‚„„f(®)‘. 

Gleichung (26): In der oberen Zeile „‚h”“ statt „Ah“. 

Gleichung (27): „R,_1,5+1 statt „R,_ 1,941 - 

Gleichung (34): ie unteren Zeile ‚(12)‘ statt 

„(1l)”. 


Gleichung (38): Man füge hinzu ‚‚(r = 2,3)‘. 
Textzeile unter (38);. „‚Iterationsformeln‘ statt ‚‚In- 
tegrationsformeln‘“. 

Textzeile unter (41): Hinter „» =1‘ füge man 
“m=0" ein. 
Gleichung (46): ‚fr‘ statt „f‘. 
Gleichung (53): Das Ende der 3. Zeile lautet 
(m+1-+2 0) “ 
4 »Y (2r) 8 
Gleichung (57): „Kg statt „ko. 
Gleichung (60): „|ß,, „|“ statt „IA, pl“. 
Gleichung (61): „A,_, statt „KRuoı- 


Gleichung (64): Der Anfang der 2. Zeile lautet 


v0 
4. Textzeile unter (64): eg |" statt [Br Be 


lebhafte Entwicklung eingesetzt hat, nicht nur in den. 


lische übersetzte 
Buch wird auch in dieser Auflage ein gern benutztes 
es mit dazu anregen, daß 
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na a ae 


ZuG.Schulz: „Zwei Hilfssätze aus der Kom- 
binatorik“. Z. angew. Math. Mech. 28 (1948) 
274—275. h | 

Herr G. Schulz hat in dieser Zeitschrift die fol- 
genden beiden Sätze bewiesen: 

I. Gegeben ist eine Strecke, die in » gleiche Teile 
(Einheiten) geteilt ist. Dann ist die Anzahl der Mög- 
lichkeiten, % nicht unterscheidbare, ! Einheiten lange 
Strecken so auf die gegebene Strecke zu verteilen, daß 
die Endpunkte jeder der k verteilten Strecken mit Ein- 
teilungspunkten zusammenfallen und keine . 
deckungen stattfinden (l ganzzahlig) 


sw N: Ki r) ah en a). 


'TI. Gegeben ist: ein Kreis, dessen Umfang in n gleiche 
Teile (Einheiten) geteilt ist. Die Einteil unkte 
seien unterscheidbar. Dann ist die Anzahl der Möglich- 
keiten, k nicht unterscheidbare, / Einheiten lan; Beh 
bögen vom gleichen Radius wie der gegebene 
auf diesen zu verteilen, daß die Endpunkte jedes der k 
verteilten Bögen mit Einteil unkten zusammen- 
fallen und keine Überdeckungen stattfinden (} ganz-. 
zahlig) durch die Formel 


u ar BET rl Re) 


gpgeben. 

Sieht man von der geometrischen Einkleidung ab, 
so handelt es sich dabei um folgendes: Gegeben sind 
n verschiedene Elemente (die Einheitsstrecken bzw. 
-Kreisbögen), für welche eine bestimmte lineare bzw. 
zyklische er a ist. Es wird die Menge 
M: derjenigen }-gliedrigen Kombinationen dieser Ele- 
mente betrachtet, deren Elemente eine Sequenz bilden, 
d.h. in der gegebenen Anordnung aufeinanderfolgen. 
Aus M sollen nun auf alle möglichen, Arten k elemente- 
fremde Kombinationen ausgewählt werden. Man hat 
es also hier mit Kombinationen von Kombinationen zu 
tun, wobei jede dieser beiden Arten von Kombinationen 
einer einschränkenden Bedingung unterworfen ist. 

Kombinatorische Anzahlbestimmungsprobleme las- 
sen häufig verschiedene äquivalente Deutungen zu und 
es kommt gelegentlich vor, daß sich darunter eine Deu- 
tung findet, durch welche das Problem auf eine der 
Grundformeln der Kombinatorik zurückgeführt: wird. 
So ist es auch beidem ProblemI. Bekanntlich ist 
die Anzahl der Permutationen von & gleichen Dingen « 
und ß gleichen Dingen b gleich 


Een a W + P) ER a (3). 


Nimmt man als Dinge a die Bedeckungsstrecken der 
Länge / (also & = k) und als Dinge b die unbedeckten 
Einheitsstrecken (also 6 = n— kl), so liefert die An- 
wendung von (3) unmittelbar die Formel (1). Es be. 
steht: eben ein ein-eindeutiges Entsprechen zwischen 
den im Satz I betrachteten Bedeckungen und den ge- 
‚nannten Permutationen, auf Grund dessen z. B: der 
Bedeckung (n = 13, k= 3, 1= 3) 


a a a a 


umkehrbar eindeutig die Permutation babbaab zuge- 
ordnet ist. 


reis 80 


liegenden Lücken, so zeigt sich, ee. ic 


Anzahl der geordneten additiven Zerlegungen 

n—klink+ 1 nie ve, ganzzahlige 

den ist 4=1+2+0+]). zu 
Zu einer expliziten Formel für das Paper N 


la; ch durch folgende 
run aut A Liegt er, so 


bleiben für die übrigen k—1 Bögen noch 17,” 


Möglichkeiten. .Man erhält also für alle k Bögen zu- 


sammen zunächst die Zahl von a - 8{"- ® Möglich- 
keiten. Da jedoch hierbei i jede derselben k mel gezählt 
ist, folgt 


x “= —H+k-1 
ze = = (Mr nn .\ 
o-M + me 
TAT —EN 


Zu demselben Wert gelangt man durch Einsetzen von 
(1) in (2). | 

Auch ie a Verallgemeinerung vonl 
2 sich auf Permutationen mit Wiederholung zurück- 
ühren: 

I*. Gegeben ist eine Strecke, die in n gleiche Teile 
(Einheiten) geteilt ist. Auf wievielfache Weise kann 
man k, Strecken der I,, k, Strecken der Länge er 
k,, Strecken der Länge l,, von 
nicht unterscheidbar sind, so auf die gegebene Strecke 
verteilen, daß die Endpunkte jeder der verteilten 
Strecken mit Einteilungs zusammenfallen und 
keine Überdeckungen stattfinden ? 

Ey: handelt sich hier um die Ta a k, 
gleichen Dingen a (Strecken der Länge l,, k, gleichen 
Dingen b (Strecken der Länge 1,) usw.undn— kl 20 
@=1,2,...,») gleichen Dingen p (unbedeckte 
Strecken der Seen 1). Die bekannte Formel liefert 
hierfür die Anzah 


m, 7, „BER 
Kla...dn Kalkal...k,ı (m kl)! 

Genau wie oben läßt sich auch die Lösung des ent- 

sprechend verallgemeinerten Problems II*, auf dessen 


Formulierung hier verzichtet werden kann, angeben, 
nämlich 
(N— N) 


{n) N 
VA  " kk—1k,...%k, 
BT 1 Kisısdy (6). 
_n. R—Zkik+ Zh—D! 
kılkt...k'im— Zkih)! 


E. Schönhardt. 


Die Zuschrift hat Herrn Prof. Dr. G. Schulz 
vorgelegen. W. 


Stuttgart. 
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erkannt werden, anderenfalls wird das folgende Quartal noch geliefert. Verantwortlich für den Anzeigenteil: CurtF,W.Schreiber 

Berlin. Druck: Treuhandbetrieb Julius Beltz ‚„ Langensalza. 57273/4022 — Nr. 994. — Veröffentlicht unter der Lizenz-Nr. 248 
der Sowjetischen Militär-Verwaltung in Deutschland. 
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